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We present skeleton-motivated evaluation of QCD observables. The approach can be applied in 
analytic versions of QCD in certain classes of renormalization schemes. We present two versions 
of analytic QCD which can be regarded as low-energy modifications of the "minimal" analytic 
QCD and which reproduce the measured value of the semihadronic r decay ratio r T . Further, we 
describe an approach of calculating the higher order analytic couplings Ak (k — 2, 3, . . .) on the 
basis of logarithmic derivatives of the analytic coupling Ai(Q 2 ). This approach can be applied in 
any version of analytic QCD. We adjust the free parameters of the afore-mentioned two analytic 
models in such a way that the skeleton-motivated evaluation reproduces the correct known values of 
\Q , rv and of the Bjorken polarized sum rule (BjPSR) db(Q 2 ) at a given point (e.g., at Q 2 = 2 GeV 2 ). 

We then evaluate the low-energy behavior of the Adler function d v (Q 2 ) and the BjPSR db(Q 2 ) in 
the afore-mentioned evaluation approach, in the three analytic versions of QCD. We compare with 
the results obtained in the "minimal" analytic QCD and with the evaluation approach of Milton et 
al. and Shirkov. 

Changes in v3: the values of parameters of analytic QCD models Ml and M2 were refined and 
the numerical results modified accordingly; the penultimate paragraph of Sec. II and the ultimate 
. paragraph of Sec. Ill are new; discussion of Figs. 4 was extended; new references were added. 
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I/-) . I. INTRODUCTION 

<N ; 

In perturbative QCD (pQCD), the coupling parameter a(Q 2 ) = a s (Q 2 )/i: [where: Q 2 = -q 2 = -(q ) 2 + q 2 ] is 
obtained on the basis of the perturbative /3-function which is a (truncated) polynomial of a. As a consequence, a(Q 2 ) 
has Landau singularities in an infrared space- like zone {Q 2 > 0), and therefore these singularities are unphysical. This 
■^5 problem was fully recognized and a solution found about ten years ago by Shirkov and Solovtsov The solution 
found was minimal in the sense that the analytization a(Q 2 ) i— > Ai(Q 2 ) was performed by removing the Landau-cut 
II ■ singularities, while keeping the singularities on the time-like axis unchanged. Further, completely analogous minimal 
Cln! analytization was performed for the higher powers a k i— > Ak (k > 2) and this replacement was performed term- 
by-term in the simple truncated perturbation series (STPS - in powers of a) of observables by Milton, Solovtsov, 
Solovtsova, and Shirkov 

US II ("Analytic Pertubation Theory" - APT). 1 The resulting series have in general better 
£> . convergence behavior and much less sensitivity under the variation of the renormalization scale (RScl) and scheme 
(RSch) . We will call the analytic QCD model based on the afore-mentioned analytic coupling the "minimal analytic" 

(MA) model [i— ► A < { AA \Q 2 )], and the afore-mentioned evaluation approach (involving the truncated analytic series) 
the APT-evaluation approach. 

The MA coupling ^ MA ^(Q 2 ) contains just one adjustable parameter, the QCD scale A. Reproduction of the 
measured values of the higher energy QCD observables (\q 2 \ > 10 GeV 2 ) fixes the scale parameter to the value 
A(„ /=5 ) ~ 0.26 GeV, corresponding to A(„ /=3 ) w 0.4 GeV. However, then the well-measured value of the massless 
strangeless semihadronic r-decay ratio r T (AS — Q,m q = 0) = 0.204 ± 0.005 @, H, [1] (cf. Appendix |E|) cannot be 
reproduced [H unless large values of the u, d and s quark masses are introduced (m q « 0.25-0.45 GeV) [l(| and 
the threshold effects become very important. One may want to avoid introduction of such large quark masses, by 
modifying the MA model at low energies while keeping the analyticity of Ai(Q 2 ) in the non-time-like region. In this 

work we introduce two somewhat different modifications A_4i(Q 2 ) (Ai = Af* 1 ^ + A.4i), both having power- like 
behaviors. We construct in a systematic way the higher order couplings Ak{Q 2 ) based on the logarithmic derivatives 
of A\(Q 2 ). Further, we construct a skeleton-expansion-motivated algorithm of evaluation of QCD observables, which 
can be applied in any analytic version of QCD and in a large class of renormalization schemes. For such an evaluation, 
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we have to know the first few coefficients of STPS and all the leading-/3o coefficients of the full perturbation series. We 
believe that the inclusion in this evaluation of the light-by-light contributions, if they contribute, should be avoided. 
Such contributions have a different topological structure and their evaluation should be performed separately in most 
evaluation (resummation) methods - see, for example, Ref. (llj . Some of the main results of the present work were 



published by us in a summarized form in Ref. 12 1 



In Sec. II, we explain the main features of the analytic versions of QCD (anQCD), we present the known MA 
model, and propose two versions of modified MA - the models 'Ml' and 'M2' [i-> _4 M1) (Q 2 ), ^ M2) (Q 2 )]. In Sec. Ill, 
we introduce the higher order couplings Ak{Q 2 ) {k > 2) in a way that can be applied in any version of anQCD, by 
imposing on them specific natural behavior under the change of scale Q 2 and of RSch. In Sec. IV we then present 
an algorithm which allows us to evaluate any QCD observable in any version of anQCD, an algorithm motivated by 
the skeleton expansion. In Sec. V we fix the free parameters in the Ml and M2 anQCD couplings Ai 1 \Q 2 ) and 
4 M2) (Q 2 )m such a way that the afore-mentioned skeleton-motivated approach gives us the measured values of r T and 
of the Bjorken polarized sum rule (BjPSR) db(Q 2 ) at Q 2 = 2 GeV 2 . We then present the resulting low-energy curves 
for the F-channel Adler function d v (Q 2 ) and of the BjPSR db(Q 2 ) in the skeleton-motivated approach, in the anQCD 
versions MA, Ml, M2. We investigate the RScl- and RSch-dependence of the numerical curves, and in the MA-case 
we compare the results of db(Q 2 ) obtained by our skeleton-motivated evaluation approach with those of the APT 
approach of Refs. [HQ. Numerical calculations were performed using Mathematica (l3j . In Sec. VI we present our 
conclusions and prospects for further work in this direction. Appendix IA1 contains details of the coefficients appearing 
in the evaluation method. In Appendix [B] we present another evaluation method that is even more closely related to 
the skeleton expansion. Appendix [Cl contains a derivation of the leading skeleton (LS) characteristic function of the 
BjPSR, and relations between the space- like and time-like formulations for the LS-term. Appendix [D] is a compilation 
of expressions of some coefficients used in this work, and Appendix [E] describes an extraction of the experimental 
value of r T (AS = 0, m q = 0). 



II. MINIMAL ANALYTIC QCD AND TWO EXTENSIONS OF IT 

The perturbative QCD coupling a(Q 2 ) = a s (Q 2 )/ir in the space-like region [Q 2 not in (— oo, 0)] has the scale 
dependence governed by the renormalization group equation (RGE) 



da{\nQ 2 ;l3 2 ,...) 
d\nQ 2 



= - ]T/V 2 a J (In Q 2 ;/? 2 ,...) , (1) 



where the first two coefficients /?o = (1/4)(11 — 2nf/3) and (3\ = (1/16)(102 — 38n//3) are scheme-independent in 
mass-independent schemes, and the other coefficients flj (j > 2) characterize the RSch. In practice, the above sum 
is truncated at a certain j max where j max — 1 is the loop level. The perturbative RGE ([I]) has a standard iterative 
solution in the form 

«(« 2 ) = EE^^. ( 2 ) 

fc=i e=o 

where L = ln(Q 2 /A 2 ) and K^i are constants depending on the coefficients and on the choice of the scale A. If the 
conventional ( "MS" ) scale A = A [3, [l5| is used, the coefficients Km are 

K w = I/fa; K 2Q = 0; K 21 = -ft/flo? 

K 3 o = -Pl/Pl + /3 2 //? 4 ; K 31 = -K 32 = -P 2 /Pl; ... (3) 

Further coefficients Km, up to k — 6, are given in Appendix [Dl The coupling a(Q 2 ), Eq. ([2]), has nonanalytic structure 

along the time-like axis Q 2 (= —q 2 ) < 0. In addition, it has singularities in the space-like region < Q 2 < A , which 
are formally the consequence of the (truncated) power expansion structure of the beta- function on the RHS of Eq. |1| . 
Application of the Cauchy theorem to function a(Q 2 ) in the Q 2 -plane gives then the following dispersion relation for 
a: 

am =ir f^vi, (4) 
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where p^ (a) is the (pQCD) discontinuity function of a along the cut axis in the Q 2 -plane: p^' t \a) — lma(—a—ie). In 
the integration, 77 is positive (77 — » +0 can be taken), reflecting the fact that the corresponding contour integration path 
avoids entirely the singularities of a(z) in the complex plane, including the singularity at z = —a = A 2 [cf. Eq. l[2"])]. 

By special relativity and causality, observables are analytic functions of the associated physical momentum squared 
q 2 = —Q 2 in the Q 2 -plane with the time- like axis {Q 2 < 0) excluded. Since QCD observables are functions of the 
invariant coupling a(Q 2 ), both should have the same analyticity properties. The singularity sector < Q 2 < A 2 in 
a(Q 2 ), Eqs. ((2]) and {4|, is therefore nonphysical. The most straightforward rectification of this problem is to eliminate 
that sector from the dispersion relation ((4]) while keeping the pQCD discontinuity function p( pt )(<7; P2, ■ ■ ■) unchanged 
on the time-like axis er > [l| , thus leading to the specific "minimal analytic" (MA) coupling 

In practice, truncated series ([2]) can be used to obtain the discontinuity function pf^ju) and thus the coupling ([5]). 
Prescription J5l) was investigated from calculational viewpoints in Refs. [r| E3, EH . There exists a practical iterative 
solution [IB, uj\ to RGE {T]) based on the Lambert function [IB] ■ This solution is an expansion of a different form than 
@ • When the number of terms in the Lambert-based expansion and in expansion ^ increases, the two solutions for 
^( MA ) conver g e to the exact numerical solution rapidly for all Q 2 , but the Lambert-based expansion converges faster. 
When fe max > 4 in @, the corresponding solution A^^^Q 2 ) differs in MS RSch from the exact numerical solution 
by less than one per cent for all Q 2 > [17]. In the present work, we will use expansion ([2]) with fc max = 5 or 6. 

Other types of analytization of a can be performed by focussing on the analyticity properties of the beta function 
[20L l2lj]. or by subtracting certain power correction terms l/(<5 2 )" from the MA coupling Ai [23]. For a review of 
various models, see Ref. [231 ] . 

In general, the discontinuity function can be different, and the analytic coupling must have the form 

■>2\ 1 f 00 dapi(a) 



Ja=0 i a + Q ) 

where pi{u) = ImAi (— cr — ie) . Relation §E§ defines an analytic coupling in the Q 2 -plane excluding the time-like 
semiaxis — s = Q 2 < 0. On this semi-axis, it is convenient to define the time-like coupling [241. 12H |26| 

2ti( S ) = f r~" Km*') ■ (7) 

2tt J_ s+%t a' 

The integration here is in the Q 2 = a 1 plane avoiding the (time-like) cut a' < 0. The relation between Ai(Q 2 ) and 
2li(s) is the same as the relation between the (vector channel) Adler function Dy(Q 2 ) and its time-like analogue, 
the e + e~ hadronic scattering cross section ratio Rv(s). Therefore, while the leading QCD correction to Dv(Q 2 ) in 
anQCD is Ai(Q 2 ) [— the anQCD analogue of a(Q 2 )], the leading QCD correction to Rv(s) is 2li(s). The following 
additional relations hold between Ai, 2li and p\ in any anQCD: 

1 f°° rlrr 

2tr( S ) = - J — Pl (a) , (8) 
"' Sli(<7) = --pi(a) . (10) 



diner 



The MA coupling ([5|) contains only one free parameter, the value of the (MS) scale A, which is not equal to the value of 
A in pQCD, but has to be adjusted so that the measured values of QCD observables be reproduced. By introducing and 
using a specific evaluation method within the MA QCD, the authors of Refs. @, H, 0] reproduced the measured values 
of the higher energy QCD observables (\q 2 \ > 10 GeV 2 ) when the scale parameter A had the value A(„ =5 ) w 0.26 GeV 

(where nf is the number of active quark flavors). This corresponds to A(„ /=3 ) ea 0.4 GeV. However, the measured 
value of the massless part of the semihadronic strangeless r-decay ratio r T (AS — 0, m q = 0) = 0.204 ± 0.005 @, H, [1] 
[cf. Appendix [El Eq. (|E6[) ] cannot be reproduced with such values of A Q unless large masses of u, d and s quarks 
are introduced (m q « 0.25-0.45 GeV) [k| and the mass threshold effects become central. 

The above consideration motivates us to introduce low-energy modifications of the MA coupling. Modifications, 
although simple, introduce additional parameters which have to be fixed by requiring reproduction of the measured 
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values of low-energy QCD observables, including of r T . One possible modification is inspired by the well measured 
@, [H isovector hadronic spectral function R v (s). At low energies (s < 1 GeV 2 ) it is dominated by the p-resonance 
(Mp = 0.776 GeV), which, in the narrow width approximation, can be represented as a delta function S(s — M 2 ) 
|27| . This is in the spirit of the Vector Meson Dominance (VMD) . If we assume that the s-dependence of the time- 
like quantity Rv(s) is at least qualitatively described by the first order time-like coupling 2li(s), Eq. ([8]), then the 
afore-mentioned delta-like structure should appear in it. This then leads to the following ansatz (model 'Ml'): 

2l^ M1) ( s ) = c f Ml6(s-Ml)+k Q(j4l-s) + e{s-M 2 )^ A) (s) , (11) 

2 2 2 2 

where c/, ko, c r = M r /A, Co = M /A arc four dimensionless parameters of the model; Q(x) is the Heaviside step 

2 

function (+1 for x > 0, zero otherwise). In this model, the MA behaviour of 2ti(s) at low energies s < M has 

been replaced by a constant (ko) plus a delta function (at s = M r < M ). The more literal application of the VMD 
approach results in ko = —1 This is so because Rv(s) = 1 + 2li(.s) + 0(2l 2 ), and Rv(s) — * when s — > 0, 

implying 2ti(s) — > — 1. However, such a model appears to restrict the low energy behavior of 2li(s) and of Ai(Q 2 ) too 
severely, especially if we want to impose the condition of merging Ai(Q 2 ) of the model with A ( ™ A) (Q 2 ) at high Q 2 . 
As a consequence, values of various unrelated low energy observables, such as Adler function (or: r T ) and Bjorken 
polarized sum rule, cannot be reproduced simultaneously in such a model. Therefore, unlike the choice ko = — 1 in 
Ref. [28|], we keep here the constant fco iu Eq. (jlip free. Applying transformation ([9]) to expression (fTTj) gives the 
space-like analytic coupling of the model 

4 M1 V) = 4 ma >(q2 ) + A 4mi) (q2) (12) 



AA ( r>m = -- 1 7" 1 +c f - x2 -i n ° 2X , (13) 

W CT=0 i^ + Q 2 ) '(Q2 + M 2 ) 2 \Q 2 + MT 



where the constant df is 



™ J All ° 



The coupling P^|) - (fT5)) can also be rewritten in a somewhat different, but equivalent, form 

A (Mi) (n2 , M 2 r Q 2 Ml Q 2 1 [°° dap^\a)(a-M 2 ) 

•^l IV ) = c f =q ' H —5 / 7 TTv?\ 

(Q 2 + MZY (Q 2 + M n ) (Q 2 + Mn) n Ja=M n cr(a + Q 2 ) 



(14) 



(15) 



In general, this coupling differs from the MA coupling ([5]) by terms A^l^ M1 ' ) ~ A 2 /Q 2 . However, we will choose to 

require SA^ 1 ^ ~ A /Q 4 , i.e., that Ml effectively merge into MA at higher energies, as we did in Ref. [28|. This 
condition eliminates one of the four new parameters, for example kg: 



C r Cf 



— =2/ W + - "Pi W 



(16) 



Since the presented version of Ml merges with MA at higher energies, the value of the scale parameter A remains 
practically unchanged, A(„ /=3 ) = 0.4 GeV, and the model contains only three dimensionless parameters c/, c r and 
Co- 

Another, somewhat simpler, modification of the MA coupling consists in adding a constant value (c„) in the low- 
energy region of the MA time-like coupling (model 'M2'): 

2ir'(s) = 2ir A) ( S ) + c„6(M 2 - S ), (17) 

Ar\Q 2 ) - A^ A) (Q 2 ) + c v M L 2 , (18) 

(Q 2 +M p ) 

2 2 

where c v and c p = M p /A are two dimensionless parameters of the model. For simplicity, we will assume that the 
scale parameter is unchanged: A( n/=3 j = 0.4 GeV. The resulting additional term cx 1/(Q 2 + M p ) in Ai(Q 2 ) can be 
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interpreted, or motivated, as the leading power-like modification (oc l/Q 2 ) of the MA coupling such that the condition 
\Ai(Q 2 = 0)| < oo is preserved. The latter condition is regarded as desirable in our approach developed in Sec. HVl 
because the so called leading-skeleton resummation of observables remains finite in such a model. 

Model Ml was motivated by simulating roughly the p-resonance contribution in the one-loop expression for Ry(s), 
via a VMD narrow width approximation ansatz in 2ti(s). However, this was only a motivation for the construction 
of an explicit form of 2li(s) as the starting point of the model, and the higher-loop contributions 2lfc(s) and Ak(Q 2 ) 
(k > 2) are then constructed on the basis of this 2ti(s) (see the next Section). The approximation of the p-resonance 
is then expected to get worse at higher loop level. Another possible approach, which we will not follow here, would 
be to refine (retroactively) 2ti(s) so that higher-loop evaluations of Rv(s) give us a given specified approximation of 
the p-resonance at low energies. A similar approach could possibly be followed also in M2. In general, reproduction 
of the correct low-energy behavior of time- like observables such as i?y(s) represents a difficult problem. In this work, 
we will follow a more modest approach - in Sec. |V]we will fix the free parameters of models Ml and M2 by requiring, 
at loop-level three or four, the reproduction of the central experimental values for the Bjorken polarized sum rule 
db(Q 2 ) at two (in Ml) or one (in M2) values of scale Q (> 1 GeV), and the reproduction of the measured value of 
ry(AS = 0). 

All the versions of anQCD presented here are infrared finite, i.e., the zero momentum limits *4i(0) = 2li(0) are 
finite. 



III. ANALYTIZATION OF HIGHER POWERS OF THE COUPLING PARAMETER 

In the previous Section, a few of the possibilities of constructing the analytic version Ai(Q 2 ) oia(Q 2 ) were presented. 
For evaluation of QCD observables, the analytic versions of higher powers a k (Q 2 ) are needed as well. For that, there 
is no unique way of constructing the correspondence a k <-> Ak- In the MA QCD, one possibility is to apply the MA 
procedure (0 to each power of a Q : 

,(MA)^ n2 ^ _ 1 f°° d(J j pt ) 
7T Jo CT + Q 2 



a k (Q 2 ) - Ar A >(Q 2 ) = - -?7rf >) (fc = l,2,...), (19) 



where pj^ — Im[a fc (— a — ie)], and a is given, e.g., by Eq. (|2|). Other choices would be, e.g. a k i— » A k ,A k 2 A 2 , 
etc. With construction (|19[) . it was shown [l6| that the RGE's governing the evolution of life's are identical to those 

governing the evolution of a k, s in pQCD when the replacements a J > Aj A ^ are made [cf. Eq. (JT])] 

a „ (MA)/ 2^ An** 
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d In n 

p 3=2 

a2 /(MA) / 2\ jm« 

2 % ] = kJ]/i,^ + t-l)fl w W = *(* + l)/«^S ) (|i , )+- , etc. (20) 

The reason for this lies in the fact that a k , and consequently p^\u), fulfill analogous RGE's. Further, the renor- 
malization scheme (RSch) dependence in pQCD, i.e., dependence of a k of (3j (j > 2), is known [29] (cf. also [30l|). 
the same dependence holds for the discontinuity functions pl pt ^ (a, /3 2 , ■ ■ ■) and thus for the MA couplings (fl9|) the 
analogous dependence via a° <-> ^ MA) is obtained (k = 1,2,...): 

dAr){p2) = ^V) + SC(^) + OiAgP) , (21) 



d(3 2 (3 k+z ^ ' 3/V 



d-4 MA V 2 ) _ k (MA) / 2\ k Pl .(MA), 2n 1 n( a(MA) 

d(3 3 - 2f3 k+s ^ ' 6f3 2 k+i ^ ) + U ^ +5 



g4^V) _ ±_ ,(MA) f 2 , Q( .(MA), , , 

The RGE-type relations (j20|) - ([23]) . valid in the MA QCD, imply the following important property: If the evaluation 
of a space-like QCD observable quantity T>(Q 2 ) is based on the analytization of STPS of that quantity according to 

the rule a fe (p 2 ) i— > A^ A \[i 2 ) (k > 1), then the evaluated value of V(Q 2 ) has a dependence on RScl /i and on RSch 
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(Pj, j > 2) which is suppressed systematically. The suppression gets stronger as the number of terms increases, just 
as in pQCD. The precision 0{A^ AA " > ) corresponds in pQCD to the precision 0{a n ). 

Having the STPS with terms up to ~ a nmax (n max = ^m), as well as its analytized analog 

4tps(Q 2 ) = «(^ 2 ;/32,...) + E^-ia n ( M 2 ;/3 2 ,...) , (24) 
^in m) (Q 2 ) = A 1 (fj, 2 ;^,...) + Y^d n - 1 A n ( f x 2 ;0 2 ,...) , (25) 

n=2 

it is then enough to include in the evolution rules (|20p and (f2"T]) - (|2"3"f (for k = 1 only) terms of up to An m on 
the RHS. Then the analytized evaluated values T> an ,(Q 2 ) will have the RScl- and RSch-independence precision 

dV ( ^\Q 2 )/dX ~ A nm+1 (X = In// 2 ,/?,-) which has its perturbative analog 8V^ s (Q 2 )/dX ~ <z n » +1 . 

In view of these considerations, we propose to maintain evolution relations (|20| (for k = 1) for any version of 
anQCD, including models Ml and M2 of the previous Section, truncating them as just mentioned: 



dA 1 ( f i 2 ;P 2 ,- 
d In /i 2 

d 2 A 1 (fi 2 ;p 2 ,- 
<9(ln/x 2 ) 2 



— —P0A2 — ■ ■ ■ — /?„ m _2^4n m , 

= 2$A 3 + 5/5 /5i^ 4 + ■ • • + n^An^ , etc. , (26) 



where we have altogether n m — 1 equations, and K n are the corresponding coefficients of the pQCD evolution equations. 
Eqs. (|26|) represent definitions of .A^'s (2 < k < n m ) via combinations of derivatives 9™^4i/9(ln/x 2 )". 

On the other hand, evolution equations (f2T1l - (|23ll (for k = 1) for the change of RSch remain of the same form, but 
with aforementioned truncation 

&*i(/i 2 ;ft,-) „ Las + ^ + .-. + ^A 

d{3 2 /3 3/3 2 ^ + + n - n - ' 

etc , (27) 

where we have altogether n m — 2 equations, and fc„ are the corresponding coefficients of the pQCD evolution 
equations. Eqs. (|27p are, in contrast to Eqs. (|26p . not definitions, but in general approximations for the evolution 
under RSch-changes. The RSch-dependence of Ai([i 2 ) is treated in more detail later in this work. 

On the basis of Eqs. (f2"6"]l - (|2~T|) , expressions for the (truncated) derivatives dAk/dX, for k > 2 (X = ln/i 2 ,/^), can 
be obtained. 

In our approach, the basic space- like quantities are Ai(/j, 2 ) of a given anQCD model (e.g., MA, Ml, M2) and its 
logarithmic derivatives 



_ (-I)"' 1 d"-Ui(/i 2 
flj _1 (n-l)! ^(ln^ 2 )"- 1 



An(^) = 7^TT 7TT 5n 77/ > (n = 1,2,3, ...) , (28) 



whose pQCD analogs are 



a n (/x 2 )= ^7 n * °L , (n = 1,2,3,...). (29) 
/?□ (» — 1)! a(ln/i- i )" 1 



The quantities (^4i(/i 2 ), ^(z^ 2 ), -43 (/^ 2 ), ••■); all derived from _4i(/i 2 ) = „4i(^ 2 ), are known functions of the space- like 
momenta /i in any chosen anQCD version in a given chosen RSch (/?2, . . .). On the basis of these quantities and the 
(truncated) evolution equations (|26p . any higher order quantity Ak{^ 2 ) (k > 2) can be constructed, in the given RSch. 
Further, (truncated) equations (f2"6"]) - (|2~T|) then give us the values of Ak(n 2 ) and of Ak{^ 2 ) (k > 1) in any other chosen 
RSch (/? 2 , /3 3 , . . .). We emphasize that in this approach, the higher order quantities Ak(^ 2 ) (k > 2) are not as basic, 
they are defined via Eqs. (|26[) for convenience of having closer notational analogy with pQCD formulas (and a k <-> Ak)- 
In these definitions (|26p . as well as in /3j-running Eqs. (|27p . we could have kept one more term (~ An m +i), in order 
to come closer to the exact analogy Ak = a k + NP for k > 2, where NP stands for nonperturbative contributions 
(nonanalytic functions of a at a — 0). 2 However, this is not necessary, as argued below. 



2 Ak = a k + NP holds exactly for the construction Eq. dl9| . i.e., the construction by Milton et al. in MA. 
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The basic analytization rule we adopt will thus be 

a„^An (n = l,2,...), (30) 

where A n and a n arc defined in Eqs. |28|) and (|29|1 . respectively. 

At loop level n max = n m , and in a chosen 'starting' RSch (/?2 , /?3, ■ ■ ■ ), the truncation ('tr') of the RGE-running of 
the pQCD coupling a(/i 2 ) is in principle via Eq. JT]) with j max = n max + 1 (a = ctt r , a n = a n ,tr)' The corresponding 
truncated A n — A n .t r are then 

A n (fi 2 ) = Sn + NP = a n (A°°) + NP + O (/3™"- 1 a"'"+™) , (n = 1, 2, . . .) , (31) 

and we assumed that we are in the class of the RSch's where fa ~ pg in the large-po limit. We recall that „4i = A\ and 
a! = a. The subscript (oo) in Eq. pip means that this is the quantity obtained by not truncating RGE beta-function 
(fTj), i.e., for j max = oo and keeping the same value of A in expansion Q as in the case of the truncated beta- function 
(i.e., jmax = "-max + 1)- The second identity in Eq. (|3ip thus shows, as an additional reference, the magnitude of error 
committed due to the truncation of the beta- function. Definitions (|26p of A n 's then imply 



An(fi 2 ) = a'V) + NP + O (/? "'"-"a"'" +1 ) (n = 2, . . . , n m ) . (32) 

n 

l (oo) 



Since the RGE-running (Q]) of a is truncated, we have a n — a?, + C(/3Q m 1 a" m+ "), and relations ([52")) remain 



unchanged when a n (y?) there is replaced by a>? oa \(fi> 2 ) ■ 

The /3j-running Eqs. (f2"T)) are also truncated, i.e., the RHS's there have errors ~ A nm +i, so that the changes of 
RSch entail additional errors. It can be verified that this effect, when going from a chosen 'starting' RSch (02,133, ■ ■ ■ ) 
to another RSch (/3 2 , /3 3 , . . .), modifies relations (pITj) to 

Ax (= A^ 2 )) = a(p 2 ) + O (/3 "»-V» +1 ) + NP , 

An (a* 2 ) = a n + (/3J»- 2 a n » +n ) + NP (n = 2, . . . , n m ) , (33) 

while relations (|32p do not get modified. We should keep in mind that there is yet another truncation involved, 
namely in the solution ||3J| of RGE fTJ) the sum over index k has in the calculational practice finite number of terms. 
In our calculations, we will take there fc max = n max + 2 (= j max + 1), which is so high that it does not affect "precision 
estimate" relations (|3"3"j) and (f3"2"| . 

For example, at loop level three (n ma x = 3), where we include in RGE ([1]) term with j max = 4 (thus fa), relations 
(J2BJ) are 



A(M 2 ) = A(M 2 ) + ^M^ 2 ) , MS) = A 3 (v 2 ) , (34) 
Po 



implying 



A 2 (n 2 ) = A 2 (v 2 ) - ^A 3 (^ 2 ) , A 3 (^ 2 ) = A 3 (n 2 ) . (35) 



The RSch (fa) dependence is obtained from the truncated Eqs. (f2"T)l and 

dA j (fj, 2 ;fa) _ 1 d 2 A J (fi 2 ;fa) ( _ 1 ~ . 2 



dfa -2(31 d(ln^ \~fa Mfl ^ 2) ) C? = l,2,...), (36) 

where Ai = Ay. These are second order approximate partial differential equations for Ai(^ 2 ;fa), A2(fi 2 ;fa), 
A 3 (/i 2 ; fa). Higher order terms (~ A4) are neglected on the right-hand side of the RSch-evolution equation [|36p. 

At loop level four (n max = 4), where we include in RGE ([TJ) term with j max = 5 (thus fa), relations analogous to 
(1351) are 



A 2 (S) = A 2 ( P 2 ) - ^A 3 (S) + ( ^ - £ ) A 4 (n 2 ) , 



A 3 (n 2 ) = A 3 (n 2 ) - \^Mf) , Mv 2 ) = A,( P 2 ) , (37) 
^ Po 
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while the changes of the RSch are governed by (approximate) relations 



d(3 2 
d(3 3 



1 



if 



5 /3i d 3 



2\(3ld{\ny?) 2 3!2/35a(ln/i 2 ) 3 
1 d 3 A 3 (p 2 ) 



312/34 9 ( ln/Lt 2) 3 



(j = l,2,...) . 



(38) 



Our approach is in a sense maximally truncating. Namely, the evolution under the changes of the RSch is truncated in 
such a way that {Q 2 ) / d(3j ~ A im +i- Further, our definition of A k 's (k > 2) viaEqs. J26j) [cf. Eqs. ([35]) and (f37|) ] 

involve short ("truncated") series which, however, still ensure the correct RScl-dependence dV^ 1 " 1 \Q 2 ) /d/i 2 ~ -4n m +i- 
Furthermore, it may seem that, for loop level three (n max = 3), the RHS of the first of Eqs. (fM|) represents only two 
perturbative terms [a 2 + (p\//3o)a 3 ] plus nonperturbative terms (NP). However, since taking j max = n max + 1 = 4 in 
RGE (fl]) as the basis for calculation of Ai(fi 2 ), it 3 is straightforward to show that the following holds: 



A 2 ( P 2 



A 2 {f) + %-A 3 (» 2 ) = aV) + V) + f aV) + 0(/3 o V) + NP 
Po Po Po 



(39) 



Completely analogous result holds at loop level 4 (n max = 4 and j max = 5) . 

In the MA QCD, in the approach of @], here Eq. (TIT))) for _4fc, a truncation is performed only in expansion @ for 

a (— » p'f t \a), apparently with fc max = n max ), and then powers of this truncated a are used to define and thus 
Ak (k > 2). This implies that in the MA QCD our Ak's (k = 2, . . .), on the one hand, and those of the approach 
of Milton, Solovtsov, Solovtsova, and Shirkov (MSSSh) [3,(1,01, on the other hand, are not the same, although they 
must gradually merge when the loop level is increased. This is illustrated in Figs. Q] and [H where the MA-coupling 



(a) , loop-level-3, k max -5, bMS 

A2, our 

A 2 , MSSSh 

10*A 3 , our 

10*A 3 , MSSSh 



Q [GeV] 




Q [GeV] 



FIG. 1: The coupling parameters A2(Q 2 ) and Ai(<2 2 ) in MA in MS RSch, with n/ = 3 and A( n/=3 j = 0.4 GeV, calculated at (a) 
loop-level=3 (and fe max = 5), and (b) loop-level=l (and fc m ax = 6). Presented are results of construction of Milton et id. (MSSSh) 



dUll, and of our construction. 



parameters A 2 {Q 2 ) and ^(Q 2 ) of both approaches are compared, for rif = 3, at loop-level (= n mM ) three and four, 
in MS and in RSch A, respectively. The Adler (A) RSch is defined later in Eqs. (J93J) [cf. Eq. ([MJ)]. For both A 2 and 
^3 one can see a decrease in the absolute difference between our and MSSSh methods when going from loop-level=3 
to 4, Fig. [1] in MS RSch, and Fig. [2] in RSch A. The decrease can be understood as coming largely from the fact that 
the perturbative part of this difference is 0(a 4 ) when loop-level=3, and 0(a 5 ) when loop-level=4. Further, inspection 



3 When the anQCD is not MA, but rather Ml or M2, RGE JT} and the (truncated) expansion J2J still remain the basis for calculation 
of the MA-part of Ai([J?), the difference between Ai(fi. 2 ) and .A^ M A ' (/i 2 ) being purely nonperturbative, cf. Eqs. 1121 , 1131 1. 1181 1. 



0.5 1 1.5 2 2.5 3 3.5 4 

Q [GeV] 



0.5 1 1.5 2 2.5 3 3.5 4 

Q [ GeV] 



FIG. 2: Same as in Fig.[TJ but now in RSch A, Eq. ((93}. 

of Figs. QJa) and[2ja) reveals that the ^-curves practically merge already at loop-level=3 if RSch is MS, but less so 
if RSch is A. An indication towards understanding this resides in the fact that the coefficient at a 4 of the difference 
between the two curves is proportional to (2/3o/?2 — this being in MS about one fifth of the corresponding value 
in RSch A (when n/ = 3). In Fig. [3] the coupling parameters Ai(Q 2 ) and ^(Q 2 ) of anQCD models Ml, M2 and 
MA are presented as functions of the scale Q, for specific chosen fixed parameters of the models Ml and M2 (see 

Sec. EJ and in the aforementioned specific RSch A. Note that we used fc max = n max + 2 in the calculation of pf^ via 
Eq. ([2|) in all cases, i.e., also in the MSSSh cases. In Fig. [3l loop-level=4 and fc max = 6 was taken (using our described 
approach). In Figs. [TJ [21 El the basis for calculation was the fc max -truncated series @ in the corresponding RSch. 



loop — level — 4 , ]^ msK -6, A RSch 



O . 4 




MA 
Ml 
M2 



0.5 1 1.5 2 2.5 3 3.5 4 

Q [GeV] 

FIG. 3: Same as in Fig. [2] but now Ai and A2 for various models (M1,M2 and MA) with specific model parameters (see Sec. fv) : 
c = 2.94, c r = 0.45, Cf = 1.08 for Ml; c„ = 0.1 and c p = 3.4 for M2; rif = 3 and A(„ 3 j = 0.4 GeV in all three models. The upper 
three curves are for _4i, the lower three curves are for 3 X An- All couplings are in RSch A, Eq. J93P . A2 is constructed with our approach. 

Even when already having anQCD coupling A\{Q 2 \ there is no unique way to merge analyticity requirements with 
the perturbative results at higher orders, i.e., Eq. (|32]) for Ak{Q 2 ) {k > 2). The latter relations are ensured by our 
definitions of Ak(Q 2 ) for k > 2 via relations (f2T)]) . but this is just one of the possibilities of addressing the problem. 
In MA the construction of Ai(Q 2 ) is very closely related to the perturbative solution a(Q 2 ) via dispersion relation 
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([5]). Therefore, it is very natural to keep that close analogy at higher orders, via dispersion relations (fT9|) . As a 
consequence, the RGE-type of relations (f2T))) are fulfilled in MA [16j . For a general anQCD model, this approach does 
not apply. Deviations of Ai(Q 2 ) and 2ti(s) from their MA values imply that the discontinuity function pi(a) deviates 

from its MA analog p < ^ t \a) — Ima(— a — it) at low values of a, cf. Eqs. (fTOf . pip and (fT7| . Therefore, there is no 
direct natural way of prescribing the low-er behavior of the higher order discontinuity functions pk (c) appearing in the 

dispersion relations of the type of Eq. (fT9|) for Ak, i.e., prescribing their deviations from p^ (a) — Ima fc (— a — ie) for 
k > 2. We define Ak(Q 2 ) for k > 2 by forcing them to obey the truncated RGE-type relations (j2l?)) . We emphasize that 
these relations define, in our approach, the couplings Ak(Q 2 ) for fc > 2. Thus, we indirectly define the corresponding 
discontinuity functions pk- This construction of Ak's is motivated also by the skeleton approach as discussed in 
Ref. [II]. Furthermore, as we will see later, this construction of Ak's allows us to suppress systematically the RScl- 
and RSch-dependence in the evaluated observables with the increasing order, because an RGE-type of analogy with 
pQCD is being preserved. 



IV. SKELETON-MOTIVATED EXPANSION 



Consider an observable T>(Q 2 ) depending on a single space-like physical scale Q 2 (= — q 2 ) > 0. Its perturbation 
expansion has the form 

V(Q 2 ) P t =a + d ia 2 + d 2 a 3 + • ■ • , (40) 

where a = a(p 2 ; fo, ■ ■ ■) is taken at a given RScl (p) and RSch (^2, •••)• As mentioned before, we will take the 
convention A = A, i.e., the MS QCD scale as the reference scale for p [cf. Eq. l[2|-([3|]. Further, we will work in the 
following classes of RSch: each [3k (k > 2) is a polynomial in rif of order fc; equivalently, it is a polynomial in /3 : 

k 

/3fc=E &fc ^o> fc = 2,3,... (41) 

3=0 



The MS clearly belongs to this class of schemes. In such schemes, the coefficients d n of expansion (|40|) have the 
following specific form in terms of /?o, as can be deduced from the scheme independence of observable T>{Q 2 ), e.g. by 
using relations of Ref. [29( : 

n 

di = c«/3 + c« d n = J2 , (42) 

i.e., each d n is a polynomial of order n in (3q and includes in general, in addition, a term with the negative power l//?o 
(di does not have it). In the MS scheme, the negative powers do not occur. 
We will now construct a separation of the series (|4"0")) into a sum of RScl-independent subseries 

oc 

V(Q 2 ) pt = P (1) (Q 2 ) pt + fc «^ (n) (Q 2 ) P t , (43) 

n=2 

with the following properties: (a) each dimensionless constant fc„ is RScl-independent; (b) each subseries T>^ (n > 1) 

is RScl-independent, and it is normalized so that = a n + 0(a n+1 ); (c) the subseries T>^ n \Q 2 ) vt contains all the 
leading-/3o coefficients of the following "rest" : 

±- \v(Q 2 ) pt - (g2) pt fc^iPf"- 1 ) (Q 2 ) pt j . (44) 

We will show that these conditions uniquely determine factors k n and perturbation expansions of all V^{Q 2 ). 
Further, we show in Appendix [B] that the above subseries, which always exist, would coincide with the expansions 
of the corresponding skeleton terms in the skeleton expansion of the observable if such an expansion existed in the 
considered RSch. 

We consider first the leading-/^ part of expansion (jit))) 

00 

^ 1} (Q 2 ) P t - a + E^'K^o] ■ ( 45 ) 
i=2 
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Under the change of RScl from /j? to /x 2 , using the notation L* = ln(/i 2 //x 2 ), we have by RGE (JTJ) 

oo 



a* + a^oi* + <(/?oi; + + at [ PtK + ^PiK + P2L* ) + 



(46) 



where a = a(/x 2 ) and a* = a(/x 2 ). Inserting this into expansion ([40]) we obtain the transformation rules for the 
coefficients ([42]) under the change of RScl. Specifically, for the diagonal coefficients the transformations are 



" * kk 



s=0 



* L 'k-s,k-s ' 



(47) 



where we use the notations cj£ = c^\n 2 ) and = c^(/i 2 ) (and = 1 by definition). Inserting expansion 
into expansion (|45[) we obtain 

^ 1} (Q 2 )pt = «* + «* + ^ [/? 2 C « + _ C W)] + 0(/33 a 4 ) . (48) 

This implies that the leading-/3o series (|45|) does not maintain its form under the change of RScl, since a new term 
a */3i( c iii — c ii ) a PP e ars at ~ a 3 . The RScl-" covariant" form, up to ~ a 3 , is then 



^(Q 2 )^ = a + a 2 0o(£ 



+ a° 



Po C 22 + Pl c ll 



/3 3 c (1) 

MO c 33 



(49) 



We now itcrativcly repeat the procedure: we insert expansion (|4"6")) into expansion (|49[) and, after some algebra and 
using relations (|47|) . obtain 



VY>(Q'% t = a*+ai f3 c^\ + a 3 /? 2 c£ 2 + ftc^ 



}2 „(1) 



J 1 ) 



+ 0(/3 o V) . (50) 
The new terms appearing at ~ a* here require the following restoration of the RScl-" covariance" up to order ~ a 5 : 



^ (1) (Q 2 ) P t = a + a 2 ftcft 



«2 (1) , a (1) 
Po c 22 + Pl c ll 



flSj 1 ) , 5 J 1 ) , fl J 1 ) 

Po c 33 + j "oPlC 22 +P2C n 



0(p> 5 ) • (51) 



This procedure can be continued to any required order. Expression (|51[) is now the RScl- covariant leading-/3o part of 
the full perturbation expansion (|40| . This means that it keeps its form (|5 1 1) under any change of RScl ix 2 . Variations of 
a = a([i 2 ) and of Ckk = Cfc/c(/x 2 ) under the RScl variation are governed by the RScl-invariance of the entire observable 
T> and of its perturbation expansion (|4"0")) . as reflected by relations (f4"6"|) and (|47[) . The additional terms appearing 
in expansion (|51[) . in comparison with the original leading-/^ series (|45[) , are subleading in /3o and represent effects 

beyond one loop involving diagonal coefficients . As shown in Appendix [B] the covariant leading-/?o expansion 
(|5ip is the expansion of the leading skeleton (LS) term in an assumed skeleton expansion of the observable T>. 

Now we subtract the LS expansion (|5ip from expansion (|40p. and the difference now involves only subleading- p\) 
terms 



D(Q 2 ) pt -2?«(Q 2 ) pt 



= fc 2 



a 2 + £ a"+ 2 ^ 



2) 



(fc 2 = 4 1 o ) ), 



(52) 



where the coefficients <i„ have a structure similar to that of d n 's (02]) 



4 2) = E « (n = l,2,...) 



(53) 



fc=-i 
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Coefficients c\ • are reiated to the original coefficients c\, by relations 



„(1)J2) 
L 10 c lj 


- c (1) 


- 


r (l) r (2) 
L 10 L 2j 


- c (1) 


5, 



fcycft (j=2,l,0,-l), 



(54) 



and coefficients 6fcj are those of the expansion of j3k coefficients (14"TT) in powers of (3q (including the case k = 1). 
Specifically, we have &fe,-i = (fc = 1,2, .. .). For k = 1, we have: &n = 19/4 and &io = —107/16, both numbers 
being RSch-independcnt. Now we repeat the previous construction, but now for the (canonically normalized) rest 
(l/k 2 )(D - P (1) ) of Eq. (5S]) instead of V (gDJ). Its RScl-covariant leading-/3 part then turns out to give 



k 2 V^(Q 2 ) pt = k 2 {a 2 + a 3 [foe™ 
Subtracting this from the rest (|52p . we obtain 



Plc { £ + p lC ^ +O(p 3 a 5 )} 



..(2 



V(Q 2 ) pt -V^(Q\ t -k 2 V^(Q 2 



n>l 



fc, - c (1) c (2) 



J 2 ) 



/3b 



d< 3) = P (c^ -b nC W)/ c W +k 4 /k 3 



(55) 



(56) 

(57) 
(58) 



where k^jkj, is a number ~ /3q which will be given explicitly below. The (RScl-covariant) leading-/3o part 2?^ 3 - ) of the 
canonically normalized expression {\/k^){J) — T>^' — k 2 V^) gives 



fc 3 2? (3) (Q 2 ) P t = ^a 3 + a 4 [^c^] +0(/3 2 a 5 )} , 
r (2) r O) _ / (2) _ , (2)x 

C 10 C ll — l C 21 °ll c ll J ■ 

Defining 

P (4) (Q 2 ) P t=a 4 + O(/3 a 5 ) , 
and following the procedure pattern, we subtract expression (|59p from expression (|56|) and obtain 
V{Q% t = 2? (1) (Q 2 ) pt + fc 2 P (2) (Q 2 ) P t + ^ 3 P (3) (Q 2 ) P t + k 4 V^ (Q\ t + Oiffia 5 ) 



(59) 
(60) 

(61) 

(62) 



where perturbation expansions for T>v> 's are given by (f5Tj) . ([55]) . (|59")) , (j6~lj) ; coefficients fc 2 and k% are given by Eqs. (j52"|) 
and [f5"7]); coefficients , c- 2 ', c- 3 ^ are given by Eqs. l|4"2"]). ([54]) . (|6"0"|) ; and an explicit expression for the coefficient £4 
is 



fc4 = c 



-,(2) 
"20 



(2) 



J2) 

C l -1 / (2) , (2)x . 1 (2) 

c< 2) Po _ 



(63) 



It is straightforward to check that all the coefficients /c 2 , k$, k 4 are RScl-independent [as are the subseries T>^'{Q 2 )]. 
Thus, identity (|62|) . obtained by our construction, represents identity (|43| to order n — 4. This construction can be 
continued to any order. 

In practice, we know only all the leading- j3o parts of the coefficients dj of observable T>(Q 2 ) Eq. ((40]) . i.e., all the 
coefficients cf)', and in addition, we usually know only one, two or three full coefficients (d\, d 2: and possibly d^). 
This implies that the first term T>^> on the RHS of identity (j6"2"]) is known to all orders, while the other terms (T>^ , 
2? (3) , and possibly £> (4) ) are known only in their truncated version. This means that the rest term in Eq. (|62p is, in 
such a case, £>(/^a 5 ), not C(/3ga 5 ). 

The perturbation expansion T>^ of the "leading-skeleton" (LS) term can be written in a resummed form [3l|, [HJ 



D (1) (Q 2 



F^)a(te c Q 2 ) , 



(64) 
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where F^(t) is the LS-characteristic function 4 which often can be written in a closed explicit form [3lj , In principle, 

F£(t) can be obtained for any space- like observable whose leading-/3o parts (c^) of all coefficients arc known. The 
value of C in (|64|) depends on the value of the reference scale A used in the RGE-running; in our convention, as 
mentioned before, we use A = A which corresponds to C = C = —5/3. 

At this point, we will turn to the question of the RSch-dependence of the (truncated) perturbation series (|6"2")l . The 
RSch independence of the series (|40p implies specific transformation rules of the expansion coefficients dj under the 
change of /3/s (j > 2) [H| 

di = dx , d 2 =d 2 - -5-(/?2 - /9 2 ) i 

Po 

d 3 = d 3 - 2d 1 ^{(3 2 - p 2 ) - -L(/? 3 - 3 ) , . . . (65) 
Po ^Po 

where the bars denote the values with MS RSch parameters fik = bk = X^fcjPo' ana - unchanged RScl. This implies, 
in view of relations (f42|) . (|54|) . ((60|) . specific transformation rules for erf,] coefficients. We will consider that the first 
term in skeleton-motivated expansion |62|) has a known characteristic function, cf . Eq. (|64|) , and that at most the first 
three nonleading coefficients of the perturbation expansion (I40[) of observable T> are known: di, d 2 , and d 3 - in MS 
RSch and at RScl /z 2 = Q 2 . Since each term in expansion (J62J) is RScl-independent, we can re-expand each V^(Q 2 ) pt 
{j > 2) in powers of a(Q 2 ), i.e., at different chosen RScl's Qj, in a chosen common RSch (p2,p3, ■ ■ ■)• The resulting 
subseries, however, will now be truncated since dj's for j > 4 are not known. This leads to the following form of the 
skeleton- motivated expansion (f6"!?|) : 

2?(Q 2 ) P t - P(Q 2 ) (TP s)+0(Po 3 a 5 ) , (66) 

2>(Q 2 )(tps) = vV{Q*)+tf\\Ql)+^\\Q>) + j^t«\\Q 2 ) . (67) 

where the coefficients t+ depend on the scale ratios Q 2 /Q 2 and the RSch parameters [3k (HT1) . and are written 

explicitly in Appendix \M in terms of the coefficients c|j' > , the latter comprising via Eq. ([H?)) the coefficients d„ of the 

original perturbation series (|40|) in MS RSch and at the RScl /i 2 = Q 2 . 

We now turn to the question of analytization of the perturbation series (|67[) . within a given anQCD model with 
known analytic couplings Ak, Eqs. (J6]), {28j-((37J). For the first (LS) term, the natural analytization procedure is to 
replace the perturbative coupling a(te c Q 2 ) by its anQCD counterpart Ai(te c Q 2 ): 5 

vW(Q 2 ) an ^V^(Q 2 )= f°° ^ F%(t) Ai(te c Q 2 ) . (68) 

.In t 



In contrast to expression (|64|) which is an ill-defined integral due to the Landau singularities of a, expression (|68|) 
is a well-defined integral in any given anQCD [unless A\(Q 2 ) diverges too strongly when Q 2 — ► 0]. We can adopt 
the viewpoint that any anQCD model is defined: (a) by a specific expression for A\{Q 2 ), and (b) by prescription 
(f6"8")) for calculation of the LS-terms of any space- like observable. The analytization of the other terms in (|6~T|) . after 
the choice of an anQCD model, i.e., of Ai(Q 2 ), can be performed in different ways. For example, the replacements 
a k (Q 2 ) i— > Ai(Q 2 ), Ai~ 2 (Q 2 )A 2 (Q 2 ), . . . ,Ak(Q 2 ) all appear equally natural at first, since the perturbative parts of 
these expressions are all the same to the order considered - cf. relations (|32p and (|3"3"|) . However, construction of the 



higher order couplings Ak (k > 2) on the basis of the anQCD coupling Ai, as presented in Sec. IIII1 suggests that it 
is the replacement 

r a k {Ql) » MQj) ^} a k {Q 2 )^A k {Q 2 ) (fe>l) (69) 



4 The superscript £ means "Euclidean", since the scales involved (Q 2 , te c Q 2 ) are space-like. 

5 A different approach to considering the perturbative LS term l|64| l was developed by the authors of Ref. [33l |. They present a novel 
version of the leading-/3o renormalon calculus, and consider that an OPE-term exists whose Q 2 -dependence is the same as that of the 
renormalon ambiguity of the perturbative LS term and that the ambiguity cancels in the sum ("PT+NP"). This sum can be presented in 
the LS form 1164 1 with the perturbative coupling a{te c Q 2 ) there replaced by a modified (but nonanalytic) coupling with one parameter. 
Since they work in the OPE framework, the latter parameter is observable-dependent. 
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that appears to be the most natural from the point of view of the requirement of the RScl- and RSch-invariance 
of the observables. Namely, Ak(ii?;(3 2 , ■ ■ -)' s fulfill, to the order considered, the same evolution equations under the 
changes of the RScl and of RSch as a k (p, 2 ;(3 2 , ■ ■ -)' s when the replacements (j69j) are performed everywhere. Further, 
the LS-analytization (|68|) of the first term T>^ of (|67|) is also equivalent to the term-by-term analytization (|69|) of the 

perturbation expansion of T>^ t , as is explicitly shown in Appendix [Bl The analytization ([55)1 of the TPS ([57)1 , which 
results in the "truncated analytic series" (TAS) 

V(Q 2 ) = V(Q 2 ) {TAS) +O((3 3 A 5 ) , (70) 

3 4 

2?(Q 2 ) ( tas) = V^ s \Q 2 )+t^A 2 (Q 2 2 ) + Y / t { i ) A 3 (Q 2 J ) + Y / tA ) MQ 2 ] ) , (71) 

3=2 3=2 

has, as a consequence, the suppression of the RScl- and RSch-dependence just as is known for the corresponding TPS 
in pQCD, but with a k i— » Ak- 



dV{Q 



2\ 



g^i = 0(p^A 5 ) (3 = 2,3,4) , (72) 



9P(Q 2 )(tas) 
d(3 k 



< 0{$t As) (k = 2,3) . (73) 



We are allowed, in principle, to vary in the TAS series (fTTj) three different RScl's Qj and 3 + 4 RSch parameters b 2 j 
and b$j appearing in (3 2 and (3 3 . One may want to have, for given chosen RScl's Qj , such a RSch that effectively only 

(2) 

the first coefficient t 2 in the beyond-the-LS contribution is nonzero. This implies various conditions involving the 
other five t^'s [Eqs. (lA3-(TA8ll]: 



tf = 4 3) =0; £4 J (T-i) 

3=2 



=> V{Q 2 ) = V^\Q 2 ) + t^A 2 {Ql) + 0([3lA b ) . (75) 

Specifically, if we choose for all three T>^(Q 2 ; p? = Q^)(tas) U — 2, 3, 4) the same RScl 

Ql = Ql = Ql = Q 2 cMC) , (76) 

the corresponding (3^ = bkjfll (fe = 2, 3) have the following Sbkj = bkj — bkj- 

5b 22 = c[l\c^+2C) , (77) 

5b 21 = c^c^ , Sb 20 = , (78) 

|^33 = cftMa + S^o^ii + C) - 5b 22 3(c$ + C) , (79) 

~«32 = cfo 4? + C^ioMo + 2fciic 1 1 ) ) - ^ 22 2c^ 1 ) - 6b 21 3(c$ + C) , (80) 

hibax = c[oC 2 o + C2b 10 c[o - Sb 21 2c[l } - 6b 20 3(c$ + C) , (81) 

-Sb 20 2c{ 1 \= 0) . (82) 



^5b 30 



-W - „(%„2 



c-j = <3 2 ;MS). Results ([7T ]> - <[82 |> are obtained by using explicit expressions (fA4|) - (fA8|) obtained 
in Appendix [2] applying to them conditions (f74"|) for the RScl choice ([751) . Specifically, result (|77|) is obtained by 
the requirement = 0; results ([75]) from the requirement tg 3 -* = 0, being zero both the coefficient at /3§ and at 

l//3o, respectively; results ([79"l) -(f82 |) are obtained from requirement 52*4 = 0, being zero all the coefficients at the 
/3o-powers (3%, 0q, 0$, l/j3 , respectively. 

Our evaluation method (fTTj) , with the choice of the scheme described above [Eqs. (|7S"j) . (|75)) . (f77j) - (j52"|) ]. emphasizes in 
the beyond-the-LS parts the role of the analytic couplings Ak(fi ) (fc > 2) constructed in Sec. Illll from the couplings 
Ai(m 2 )> Eq. (p8|) [see Eqs. ([37|)]. The couplings Ak{(i> 2 ) [k > 2) were constructed in such a way as to have, at 
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perturbative level, their equivalence with a n {fi 2 ). However, the construction in Sec. IIIII strongly suggests that the 
couplings A n {^ 2 ) {n > 2) are more basic since they are constructed as derivatives of A\{^ 2 ) which is the basic quantity 
in any anQCD model. Further, the skeleton-expansion arguments presented in Appendix [Bl show that A n {/J. 2 ) are the 
basic elements for the expansion of each term in the skeleton expansion. Therefore, a more natural choice for RSch 
(02, f3 3 ) in the evaluation method (f7T|) . with RScl's (|76|) . would be such that the resulting TAS expression is 

V(Q 2 ) = V(Q 2 ) {TAS) + 0($A 5 ) , (83) 
2>(Q 2 )(TAS) = V^{Q 2 )+t 2 A 2 {Q 2 e c ) . (84) 

To obtain the f3k'$ (k = 2,3) necessary for this result, we first re-express all Ak's (k > 2) in TAS (fTTj) in terms of 
A n 's, Eqs. (f3"T)) . Keeping the RScl's according to (f75|) . this implies that, in a general RSch {{fa, (is) expression (fTTj) 
can be re-expressed as 

V{Q 2 ) (T as) = ^ LS \Q 2 )+t 2 A 2 {Q 2 e c )+t 3 A 3 {Q 2 e c )+UA 4 {Q 2 e c ) , (85) 

where the coefficients U are certain combinations of ti k \ and are written explicitly in Appendix [2] Eqs. (|A16[1 - (| A21|) . 
Requiring the form (|84|) . i.e., 



h = U = , (86) 
implies, by Eqs. (|Al6j) - ()A21l) . that the corresponding flf. = 6jy/3g {k = 2, 3) have the following Sbkj = bkj — bkj'- 

Sb 22 = 4 ) (cf 1 ) +2C) , (87) 

Sb 2 i = c^Mo - b n cio , Sb 20 = -b w c[]j , (88) 

\sb 33 = cftc^ + 3#JC(j$ + C) - Sb 22 3{c$ + C) , (89) 

1 „ Ul)U2) _ 5 _(1)_(2) _„(1)T , ,^-(l),-(2) _ , x , „ (_«J?.),K\ _n orf 



-32 = Sio ; 4i " 2 6 " 5 io cii - cio J &22 + 3CZ& (chY - 6n) + ^22 I -3Z$ + -611 ) - *6 2 i3(c^ + C) , (90) 

1„ -(l)-(2) 5 _(1)_(2) 5 _(1)_(2) L _(i)/5 2 T », , 

-5&3i = cVc 20 ' - 2 6 io c io c ii - 2 11 i° i° + w { 2 11 ~ ~ 

+^22^10 + *&2i (-32$ + |& u ) - 5b w 3{c$ + C) , (91) 
5*630 = -§6io^oMo + S^iofencio 3 - ^ocio 3 + 5610*62! + *6 20 (-32$ + |fen^) . (92) 

In these expressions, b 2 j are the coefficients b 2 j in MS: 6 22 = 325/96, 621 = 243/32, 620 = —37117/1536 (and 
6n = 19/4, 610 = —107/16). We will apply, as a rule, our evaluation approach in the RSch (f87|) - (f92|) . i.e, where the 
resulting formula is (|8"3")) - (f84")) , and will use the RScl's ([7rJ)l with C = C = —5/3. The RSch evidently depends on the 
observable. Our starting point will be this RSch for the massless Adler function T>{Q 2 ) = d v {Q 2 ), where the STPS is 
known to a large degree of accuracy up to ~ a 4 (up to ~ a 3 it's known exactly) - we will call this RSch A ('A' for 
Adler). 6 If an observable is known in STPS only up to ~ a 3 , only formulas (|87 p - (|55|) are to be applied, as £4 is not 
known; in that case, in Eq. ([55]) the unknown rest term is O{0qA^). For example, Bjorken polarized sum rule db{Q 2 ) 
is such an observable. 

In Appendix[B] a different method of evaluation is presented, which would be an evaluation of the skeleton expansion 
itself if such an expansion existed in the considered RSch. The RSch-dependence of that method is numerically 
stronger, which may be a reflection of the fact that this expansion, if it exists, is valid only in a specific ('skeleton') 
RSch that is hitherto unknown [H, HH . 



The difference between this RSch A and the RSch A' 177t - d82l for the Adler function is small. For example, for rij = 3, the values 
are f3 { 2 A) = -18.92, /3^ A '' = -18.59; /3^ A) = -33.84, /3^ A '' = -32.72. In Ref. [TJ, we used RSch A' f77l - (gg) [with RScl's J76} with 
C = C = —5/3], and denoted there this approach as V2'. 
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V. NUMERICAL RESULTS 

In this Section, we take the position that the anQCD models Ml and M2, introduced in Sec. [TTl the form of 
A\(Q 2 ) there, Eqs. (fTS")) and (|18p . is achieved in the aforementioned "optimal" RSch (|57|) - (|92p for the massless Adler 
function d v (Q 2 ) - RSch A. We must keep in mind that models Ml and M2 change the form of A\(Q 2 ) when the RSch 
(02,03, ■ ■ ■) is changed. 7 

We will calculate numerically various low-energy QCD observables in the anQCD models MA, Ml and M2, with 
rif — 3, by using the skeleton- motivated evaluation method presented in the previous Section, Eq. (|85|). One such 
quantity is the massless Adler function d v (Q 2 ) whose pQCD expansion coefficients d\ and d 2 (in MS RSch and at 
RScl /j, 2 = Q 2 ) are known exactly [III |37|, and d^ has been estimated as a polynomial in n/ to a high degree of 
accuracy [38] (see Appendix iDl for explicit expressions of d\, d 2 , d%). The normalization of d v is taken according to 
Eq. (|4TJ|) when rif = 3. The additional light-by-light contributions [13] do not contribute when rif = 3. Further, 
the LS characteristic function F„(t) for d v (Q 2 ) was obtained in [3l[, and is given in Appendix [Cl in Eqs. (|C6[) - (|C7|) . 
Evaluation method (|85p can thus be applied by including terms ~ A4 in the case of the massless Adler function (for 
a different approach to evaluating Adler function, see Ref. [39|). The optimal RSch for the massless Adler function 
d v (Q 2 ) is then obtained by requiring disappearance of ~ ^3 and ~ A4 terms, Eq. (f86|) . where we choose RScl according 
to (176]) with C = C = —5/3. We call this RSch Adler (A), and it can be obtained from MS RSch by applying relations 
(ETJ-dlS), resulting in 

( 2 A) = -23.6074 - 16.0248/?o + 8.04784/^ , 

( 3 A} = 127.38 - 35.8577/3 - 12.8734/$ - 1.34926/$ . (93) 

The values for rif = 3 are 2 = -18.9211 and 3 = -33.8404 (in MS RSch, at n f = 3, the values are 10.0599 and 
47.2281, respectively). In RSch A, the evaluated massless d v (Q 2 ) is thus 

d„(Q 2 ) TAS = jH j F v £ (t) A^Q 2 ;^ \0 { 3 A) ) + ^M^Q 2 ) , (94) 

and the difference between the (massless) true d v (Q 2 ) and (4(Q 2 )tas is formally O(0qAs). 

The (V + A-channel) semihadronic r decay rate ratio r T is one of the best measured low-energy QCD quantities, 
its massless part for non-strange hadron production has the value r T (AS = 0, m q = 0) = 0.204 ± 0.005 [1, Q 
(cf. Appendix [El Eq. (|E6|I ). The heavy quarks (c and b) do not contribute, since r T is a Minkowskian observable, and 
the t particle cannot decay to charmed mesons because their masses are larger than m T . s Our evaluation approach 
for r T (AS — Q,m q — 0) uses the aforementioned evaluation ([M]) of the (massless) Adler function d v (Q 2 ) which is 
then inserted in the contour integral (|C8|) . The LS-part can then be written in the form (|C9|) with the time-like LS 
characteristic function (|CJ10[) - (ICJ1 1|) . The beyond-the-LS (bLS) contribution is the contour integral 

r T (A5 = 0,m 9 = 0)( bLS >=i dcf> (1 + e^f(l - e^)^-A 2 (eV T e^) . (95) 

27T J_ 7I 12 

Yet another low-energy QCD observable that we will consider is the Bjorken polarized sum rule (BjPSR) db(Q 2 ). 
Its LS-characteristic function is obtained in Appendix [Cl on the basis of the known leading-/?o coefficients [4fJ using 
the technique of [3lJ. The full perturbation coefficients d\ and d 2 for the massless db(Q 2 ), in MS RSch and at RScl 
/i 2 = Q 2 , were obtained in Rcfs. 41] (see Appendix |Pl for explicit expressions for d\ and d 2 ). For the coefficient d%, 
only the leading-n/ part (<x n^) is known exactly [401 ] : based on this, estimates of d% as a polynomial in 0q were 
performed in Ref. [12] using naive nonabelianization (NNA) i-> — 60o 43]. For the evaluation of (the massless part 
of) db(Q 2 ) we will not use estimates of the full d%, i.e., we will use method (|85|) with terms up to ^3^3 included, in 
any chosen RSch and with RScl's (|76[) with C = C = —5/3. The formal difference between the evaluated and the true 
value is then O(0 2 Aa). The experimental values of db(Q 2 ) at low Q 2 are much less precise than those of r T (AS = 0). 
At Q 2 = 2 and 1 GeV 2 they are 4(2GeV 2 ) = 0.16 ± 0.11 and 4(lGeV 2 ) = 0.17 ± 0.07 Q (for an application, 



7 When /3j's (j > 2) change, the change of A.i(Q 2 ) in general cannot be described just by running of the parameters of the model with 
/3j's, since new terms appear that depend on those parameters. 

8 The contributions of heavy quarks in Euclidean observables T>(Q 2 ), such as the Adler function, can be more important, even though 
Q 2 < m 2 — see the discussion later in this Section. 
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TABLE I: Results of evaluation of the semihadronic tau decay ratio r r (AS = 0,m q = 0) and of BjPSR d b (Q 2 — 2 GeV 2 ), in 
various anQCD models, using evaluation method (|85[) in RSch A (|93[) . The basis for calculation of P\ {&) is expansion ((2} 
at loop-level=4 (i.e., when included) and with fc max = 6. In parentheses are the results at loop-level=3 and fc max = 5 
(in that case, the d3-term of the Adler function is not included). Presented are the results of the full evaluation (leading- 
skeleton and beyond: LS+bLS), Eq. (|85[) . an d for r T (AS = 0,m q — 0) also the results of LS. The experimental values are 
r T {AS = 0,m, = 0) = 0.204 ± 0.005, d b (Q 2 = 2 GeV 2 ) = 0.16 ± 0.11 and d b {Q 2 = 1 GeV 2 ) = 0.17 ± 0.07. See the text for 
further details. 





r T (AS = 0,m q = 0) 


rv (AS = 0,m q = 0) [LS] 


d b {Q 2 = 2 GeV 2 ) 


d b (Q 2 = 1 GeV 2 ) 


MA 


0.141 (0.142) 


0.139 (0.141) 


0.137 (0.138) 


0.155 (0.155) 


Ml 


0.204 (0.205) 


0.197 (0.198) 


0.160 (0.161) 


0.170 (0.171) 


M2 


0.204 (0.206) 


0.203 (0.204) 


0.189 (0.190) 


0.219 (0.220) 



cf. Ref. [ill). The contributions of massive quarks (m C) m&) are \Sdb(Q 2 ;m q ^ 0)| < 10~ 3 for Q 2 < 2 GeV 2 |46j ]. 
thus negligible. We recall that both d v and dh are massless observables which are normalized here according to the 
convention (|40|) for n/ = 3. Although the uncertainty of the measured values of db(Q 2 ) is significantly lower at 
Q 2 = 1 GeV 2 than at Q 2 = 2 GeV 2 , we will use both central values. We expect the theoretical predictions of our 
evaluations in general to be more reliable at higher momenta Q 2 > 1 GeV 2 . 

Now we will fix the parameters of models Ml and M2. Model Ml (fl"T j) - (fT6)) has three independent parameters 
Cf, cy, Co (and A = 0.4 GeV as in MA). Requiring the reproduction of the aforementioned experimental central values 
r T (AS 1 = 0, m q = 0) = 0.204, (4(2GeV 2 ) = 0.16 and c4(lGeV 2 ) = 0.17, we obtain a solution for the three parameters, 
with the following values: c/ = 1.08, c r = 0.45, Co = 2.94. We will use these parameter values in Ml (in RSch A). 
In general, the predicted values of observables do not change a lot when cq is varied in the regime ~ 1; they change 
more when c r and/or Cf are varied. The experimental values of various higher-energy QCD observables T>{Q 2 ),lZ{s) 
(Q 2 y s ^ 10 GeV 2 ) should be well reproduced in Ml, because condition (fl"6]) ensures that Ml and MA merge at higher 
energies Q 2 ,s A , and it has been demonstrated that MA with A(„ =3 ) = 0.4 GeV (=> A(„ _ 5 j = 0.26 GeV) 

reproduces well those values Q. We note that model MA (with A = 0.4 GeV) predicts r T (AS = 0,m q = 0) « 0.14, 
which is significantly too low. 

Model M2 (| 1 T[) - (TTH|) has two free parameters c v and c p , both assumed to be ~ 1. Requiring reproduction of the 
central value of r r (AS = 0,m q = 0) = 0.204, and requiring |c p |, |c„| > 0.1, it turns out that the model then always 
predicts values (4(2GeV 2 ) > 0.19. Requiring the minimal possible value <4(2GeV 2 ) « 0.19 gives us the parameter 
values Cy =0.1 and c p = 3.4. We will use these parameter values in M2 (in RSch A). 

In Table 1 we present results of calculations of r r (AS = 0, m q = 0) and db{Q 2 = 2GeV 2 ) with our evaluation 
method (|85|) , in the aforementioned RSch A ([93|) - ([94)l and at loop-level=4 and 3, in various anQCD models: Ml, 
M2, and MA. When loop-level=4 (and fc max = 6), we used in the calculation of r T (AS = 0,m q — 0) the estimated 
N 3 LO perturbation coefficient d% of Ref. [38] for the Adler function (cf. Appendix [P|) . as mentioned earlier. In the 
case of db(Q 2 = 2GeV 2 ), when loop-level=3 or 4, evaluation formula (|85|) was used in RSch A by inclusion of terms 
up to A3 only, as the N 3 LO coefficient d% is not known there. We note that MA (with A( n/=3 ) = 0.4 GeV), with 
light quark masses m u ,m,},m s <C A (m u ,md,m s w 0), does not reproduce the well-measured experimental value 
r T (AS = 0, Tn q = 0) = 0.204 ±0.005, as already mentioned in the Introduction. This fact led us to suggest alternative 
versions of anQCD (e.g., Ml, M2). 

Now that the parameters of the presented anQCD models have been fixed, we can present various results of these 
models, evaluated with the method |85|) . In Fig. 0Ja) we present curves for the massless Adler function d v (Q 2 ) 
(with nf = 3) as functions of energy Q, in models Ml, M2, and MA. The RSch used is RSch A f93"] ) -(|93] ) . Loop- 
level is four, i.e., we include the value (3^ in our calculation for pf* , with fc max = 6 [cf. Eq. (J5J)], and use the 
estimated N 3 LO perturbation coefficient d$ of Ref. [38| (cf. Appendix |D|) . The light-by-li ght contributions, which 



11]), appear for the first 
371 ]. This sum is zero in 



have a different topology of diagrams and should probably be resummed separately (cf. Ref. 
time at ~ a 3 and are proportional to the square of the sum of the quark charges (Y^ Q f) 2 
the case nj = 3 considered here. Fig.[4jb) represents the results for the full Adler function, i.e., the V-channel heavy 
quark corrections 5d v (Q 2 ;m c ,mb) have been added there. For the calculation of the latter, we follow the procedure 
of including the ^-contributions (note that d v (Q 2 ) = (1/2)D(Q 2 ) - 1, where D is defined in [47}). The first 
seven coefficients of the low-momentum Taylor expansion for the heavy quark a 2 -contributions are calculated in [48| . 
Through a conformal mapping together with Pade improvement, as proposed in [49j . an approximant is obtained. 
The approximant reproduces the low-momentum behavior and fits very well the large-momentum expansion [Hp} for 
this quantity up to energies Q 2 w 16m 2 (see also Fig. 4 of Ref. [47j). Thus, this method can be safely used for 
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FIG. 4: Adler function as predicted by pQCD, and by our approach in several analytic QCD models (see the text): (a) the massless part 
(rif =3); (b) the full quantity, with the contribution of massive quarks included. 



the q = c, b quarks in the energy range we are interested in. 9 In the heavy quark contributions, we simply replaced 
a(Q 2 ) and a 2 (Q 2 ) by Ai(Q 2 ) and A 2 (Q 2 ) (using A = A = 0.4 GeV). The indicated ± uncertainties in the full Adler 
function curves are those c quark contributions which are proportional to A 2 - In Figs.[4ja),(b) we included the STPS's 
[truncated forms of Eq. ([4*0)) ] in MS RSch and with RScl fj 2 — Q 2 . In Fig. |4jb) we included experimental values, 
for comparison. The experimental values of d v (Q 2 ) are taken from Ref. |47| where the integral expression for d v (Q 2 ) 
in terms of the e + e _ QCD ratio R e + e -(s) is evaluated. All the values of R e + e -Js) are needed - from the two-pion 
threshold to infinity. The evaluation is based on the data compilation of Ref. [54]]. The pQCD result for R e + e -(s) 
is used in the integral where it can be trusted, and data in the rest of the energy interval. Resonances are included 
separately. In Fig. 0Jb) we can see that various anQCD models predict at low energies (Q < 1.2 GeV) values which 
are significantly closer to the experimental values than STPS's. Further, STPS's lose any predictability at Q < 1.2 
GeV, mainly because of the vicinity of the unphysical Landau pole in the pQCD coupling a(Q 2 ). 




1 1.5 2 2.5 3 3.5 4 1 1.5 2 2.5 3 3.5 4 

Q [GeV] Q [GeV] 

FIG. 5: Bjorken polarized sum rule (BjPSR) db{Q 2 ) in (a) model Ml, and (b) comparison of Ml and MA; in various RSch's and 
at various RScl's. The vertical lines in (a) represent experimental data, with errorbars in general covering the entire depicted 
range of values. 

In Fig. [3a), we present results of calculation of BjPSR db(Q 2 ) at low energies in model Ml, at loop-level=3 (and 



Some contributions from heavy quarks are not considered here as we base our analysis on the expressions of Ref. *4""ll . The relevant 
diagrams are shown in Fig. 2 of Ref. |48ll ; the contributions with internal heavy and external light quarks are not included. These 
type of (a 2 -)contributions have been obtained for the R e + e -(s) function in Refs. [5ll. [52f, |53(. We checked that these contributions, 
when translated into the corresponding contributions for d v (Q 2 ) via the usual integral transformation relating R and d v , result in 
a 2 -contributions which are an order of magnitude smaller than the heavy quark a 2 -contributions included in our curves. 
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k max = 5), in two different RSch's: RSch A J93|), and RSch B which is the "optimal" RSch for d b (Q 2 ), i.e., /3f> is 
obtained from the requirement t 3 — for d b , Eqs. (f57|) - (|85|) 

4 B) = -30.2949 - 10.4415/3 + 7.44582/?^ . (96) 

At n/ = 3 we have j3\ (rt/ = 3) = —16.0938. The analytic couplings in RSch B are obtained from those in RSch A by 
applying the looplevel=3 RSch-evolution equations ([3"6"| . In addition, we present in Fig.^a) results when the RScl in 
the beyond-the-LS terms (Q 2 ,Q 2 ) is increased from Q 2 exp(— 5/3) to Q 2 (note that coefficients t-i and £3 then change 
accordingly). We see that at low energies Q < 2 GeV the results in Ml change moderately but not insignificantly 
under the variation of RSch and RScl. For comparison, we included the curve obtained from the skeleton evaluation 
(|B21j) in RSch A [with Q\ = Q\ — Q 2 exp(— 5/3)], assuming that the skeleton expansion exists in RSch A (which is 
probably not true). We include the present experimental data, with the crosses representing the central values; the 
crrorbars extend in general over the entire depicted range of values, most of the experimental uncertainties are of 
the order of ±0.1. The experimental data were deduced from Fig. 2 of Ref. p5j . with the neutron decay parameter 
value 1*7^4 1 = 0.21158 ± 0.00048 taken from (5f|. The present experimental errors are too high to discriminate between 
various evaluation methods. In Fig. [5Kb) we compare the results for of MA and Ml. The RSch- and RScl-dependence 
of MA-results remains very weak in all the shown region. 




vertical lines in (a) represent the experimental data. 

In Fig. [6](a) we present the same type of curves for M2 model. We see that the RSch- and RScl-dependence in M2 
remains quite weak down to low energies. In Fig. [DJb) we compare the results of M2 and Ml models. Only the curves 
in RSch A are presented in Fig. [U(b). 

Up until now, we applied the (skeleton-motivated) method (f85|) for the evaluation of QCD observables, in various 
anQCD models for Ai(id 2 ), with the higher-order couplings Ak (k > 2) constructed by Eqs. (|28| in a certain RSch 
(usually RSch A) and equivalently the higher-order couplings Ak by Eqs. ([33]) [Eqs. ([57)) if loop-level=4] . There 
remains a question of how this method of evaluation compares with the APT evaluation approach of Milton et al. and 
Shirkov dlH. We recall that the APT approach was defined for the MA anQCD model, and it consists of using the 
available (NLO or N 2 LO) STPS of an observable (gO]) and replacing there a k (Q 2 ) i-> A k (Q 2 ) (MA) (k > 1), where the 
higher-order MA couplings Ak{Q 2 )^ MA ^ were constructed according to formula ([T9|) . In the N 2 LO STPS case [e.g., 
for db(Q 2 )], this reads 

2WQ 2 ) = A(Q 2 ) (MA) + di.4 2 (Q 2 ) (MA) + dMQ 2 ) {MA) ■ (97) 

The RSch is usually taken to be MS, but could in principle be any RSch. One of the differences between our and 
APT evaluation method here is the construction of the higher-order couplings Ak(Q 2 )( MA ^ of the model MA, where 
comparison with our construction has been made in Figs. [1] and [2] in Sec. lIIIl Another difference is that our evaluation 
method ([85| includes, in addition, the leading-/3o contributions to all orders. We compare in Figs.[7ja), (b) the results 
of our method §5\) and APT-method §7]) for BjPSR d b {Q 2 ), in MA model, in two different RSch's: A ([M]), and MS 
(relevant for (3 2 coefficient only). In Figs. [3(a) and (b) the renormalization scale was taken to be /j 2 = Q 2 exp(— 5/3) 
and [i 2 — Q 2 , respectively, in both beyond-the-LS terms of our approach (oc ^27-^3)1 and in APT approach. The 
RSch-change from RSch A to MS was performed in our approach according to the loop-level=3 evolution equations 
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1 1.5 2 2.5 3 3.5 4 1 1.5 2 2.5 3 3.5 4 

Q [GeV] Q [GeV] 

FIG. 7: BjPSR d b (Q 2 ) in the MA anQCD, with our evaluation method ((85) and that of Milton et al. and Shirkov (APT), in RSch 

A 193 > . and in MS; the RScl is chosen to be (a) fi 2 = Q 2 exp(— 5/3), (b) fi 2 = Q 2 . The vertical lines in (a) represent the experimental 
data. 



(l36|) . while in APT approach the corresponding values of (32 were inserted directly in ([2]) and thus in all s. For 
additional comparison, we included the skeleton evaluation (in RSch A), Eq. (|B21| . We see in Figs. [7] that in both 
our and APT evaluation approaches, in MA anQCD model, the RSch- and RScl-depcndcncc of db(Q 2 ) is very weak 
down to quite low energies. More detailed inspection reveals that our evaluation approach ([55]) gives for db(Q 2 ) even 
somewhat less RScl-and RSch- dependent results than APT approach. 

VI. SUMMARY AND CONCLUSIONS 

In this work we suggested various models of analytic QCD (anQCD), i.e., models for construction of the anQCD 
coupling Ai(Q 2 ) which is an analytic analog of the perturbative QCD coupling a(Q 2 ) = a s (Q 2 )/n. The main reason 

why we suggest alternatives to the minimal analytic (MA) model, i.e., to the coupling Ai~ A \Q 2 ) of Shirkov and 
Solovtsov [l[ , is that it cannot correctly reproduce simultaneously various higher-energy QCD observables on the one 
hand and the low-energy observable r T (scmihadronic r decay rate ratio) on the other hand, unless large masses of 
U, d and s quarks are introduced [|J|. The described alternative models (Ml and M2) have Ai(Q 2 ) with additional 
dimensionless parameters in it, which can be adjusted in order to modify the behavior at low Q 2 . Furthermore, we 
presented, for any anQCD model, an algorithm which allows construction of higher-order analytic couplings Ak(Q 2 ) 
(k > 2) which are the analytic analogs of a k (Q 2 ). In addition, we presented a method of evaluation of Euclidean QCD 
observables in anQCD models, a method which is (partly) motivated by the so-called skeleton expansion structure 
but does not depend on the existence of such a skeleton expansion. The evaluation method sums up all the leading-/^ 
contributions (LS: leading-skeleton) and adds those contributions beyond the LS which are known by the knowledge 
of a first few perturbation expansion coefficients of the considered observable. We tested this evaluation method, for 
three anQCD models, in the case of the Adler function, semihadronic r decay ratio, and the Bjorken polarized sum 
rule (BjPSR) at low energies. The results show in general good stability under variation of the renormalization scale 
and scheme down to low energies Q ~ 1 GeV. We further carried out comparison of our evaluation method with that 
of Milton et al. (APT) [2j,|3|, |4j, for the BjPSR, in the MA model where the latter method can be applied. The two 
methods give results which at low energies differ in general by only a few per cent for this observable. 
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APPENDIX A: RELEVANT COEFFICIENTS OF THE SKELETON-MOTIVATED EXPANSION 

In this Appendix, we present explicit formulas for the coefficients ij appearing in the skeleton-motivated expansion 
(|S7)) . which is a slightly reorganized form of expansion (fl)2")) . We consider the case when, in MS RSch (0/. = bf. = 
Yl &fe?'Po> k > 2) and at RScl [j, 2 = Q 2 , the first three coefficients in expansion PD|) are explicitly known (dj, j = 1, 2, 3), 
and all the leading- Q parts CnnPo™ of coefficients d„ (n > 1) in expansion (|42|) are known (we note that = 

in MS). The RSch (02,03, ■ ■ •) is chosen and common to all terms DW (j > 1) ; and belongs to the class of RSch's 
of Eq. ((U]). The RScl's used in the resulting truncated versions of T>^'(Q 2 ) (j > 2) are Q 2 , they may be mutually 
different as each V^(Q 2 ) (and kj) is RScl-independent. For the RSch and the RScl's we will use notations 



5b k j = b k] - b k] , (Al) 
Q 2 = Q 2 exp(C J ) . (A2) 

We then obtain for the coefficients tj of expansion (jCTT)) the following expressions, on the basis of relations (f5Tj) - (fl)3")) . 
as well as (gTJ), (g2J» and (M]): 
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Here, t[ 3) are the values of the t { f in MS RSch and with RScl [i 2 = Q 2 : 
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Coefficients cff appearing in the above formulas can be obtained directly from coefficients cj^ by using relations (|54|) 



in MS scheme (with RScl y? = Q 2 ): 
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Formulas (|A15|) . ([X9 )> - (|A14|> . lfA3 )) -([A"8 l) . with notations ([AT]) and |[A"2]) . allow us to obtain all the coefficients t> of 
the skeleton-motivated expansion (jUT)) in any chosen RSch and with chosen RScl's Q|, if we know in MS RSch at RScl 

Q 2 all the leading-/3 parts CnnPa of the expansion coefficients d n = J^o^nk^a °f observable V(Q 2 ) Eq. (JTUJ) , 



and we know exactly the full coefficients dj for j = 1, 2, 3, i.e., we know c^J for j = 1,2,3 and k = 0, . . . , j. If, on the 
other hand, we do not know ds, the above formulas can be applied for t/ for i — 2,3 only. 

When the beyond-the-LS contributions in our approach (jTTj) are expressed in terms of Ak's, Eq. (|85p. with the RScl 
choice (1751) [Cj = C] , coefficients U can be expressed in terms of the above coefficients via relations (|57|) between 
Afc's and -4 n 's. After some straightforward algebra, we obtain 
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where tj are the values of t% in MS and with RScl y 2 = Q 2 
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APPENDIX B: SKELETON EXPANSION 



In this Appendix we will construct an expression for evaluation of QCD space-like observables T>(Q 2 ) (for any 
anQCD model) which will be derived directly from the QCD skeleton expansion. Here we will take the position that 
such an expansion exists in the class of schemes with the QCD scale A 2 , = A 2 , exp(C) where Ao is the so-called l^-scheme 
scale and C is an arbitrary n /-independent constant, and with (3k of Eq. (|41|) where bkj are arbitrary constants. In 
this context, choosing the MS scale parameter C = C = —5/3 (A = A) for scaling the RScl fi 2 represents no additional 
restriction. This expansion involves in the integrands the characteristic functions F^(ti,. . ., t n ), which are considered 
n/-independent, and the (singular) pQCD coupling a(/i 2 ). We replace a(^t 2 ) by an anQCD coupling Ai(ii 2 ) in the 
skeleton integrals which makes the integrals unambiguous 
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Here, F^{t\,. . ., t n ) are the characteristic functions and have the normalizations 
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implying for the perturbative parts 

^ (K) (Q 2 ) P t = a n -[l + 0(a)] , (B4) 

where n K = 1 for k =LS, n K = 2 for k =NLS, etc. The perturbative part of Ai(/j, 2 ) is a(fi 2 ) [Ai(fi 2 ) — a{^ 2 ) + NP, 
where NP involves non-analytic in a — functions of a, cf. Eq. I|3ip ]. We will use RGE evolution series (|46f for 
expansion of a(te c Q 2 ) around a(/j, 2 ) = a 
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where T = tQ 2 e c / fi 2 , and a n are defined in Eq. (|29|) . Using expansion (|B5|) in the leading-skeleton (LS) term in (|B1[) . 
this term can be shown to have the following expansion for its perturbative part: 
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where we adhere to notations summarized in the following: 
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Requiring that the perturbative part of the LS-term absorb all the leading-/3o parts of T>(Q 2 ) pt [see Eqs. (HU|) - (H2"1) ] 
implies that 

<(-mT)«> (1) =c« (n = 0,1,2,...) . (B9) 

This, in conjunction with expansion (|B6p . implies that X>( LS ) (Q 2 ) p t is precisely D^- 1 ' (<5 2 ) P t of construction in Sec. IIV[ 
Eq. (|5Tj) , i.e., we really have for X > ' 1 ^(Q 2 ) p t the resummed form Taylor expansion of _4i(te c (3 2 ) around Q 2 is 



completely analogous to expansion (|B5|) 
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where Ak = »4fc(/i 2 ) and ^ = .4fc(/j 2 ). In the last identity we used the fact that A n 's appear on the left-hand side of 
RGE's ((26)) . which are analogous to pQCD RGE's with a„'s on the left-hand side [analogy valid up to terms 0(A nTa ) 
where n m =loop-level] when the correspondence a k <-» Ak is made. Eq. (|B10p implies for the (full analytic) LS-term 
of the skeleton expansion (|B1|) a nonpower analytic expansion 
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which is just the analyticized analog [according to the rule ((69)) ] of perturbation expansion (|B6[) and (|5l"j) . 

Now wc will investigate the beyond-the-LS contributions of the skeleton expansion (|Bip . In view of normalization 
conditions (|B4j) . it follows immediately that 



(B12) 



which is just the coefhcient k% (|52p in the approach of Sec. IIVI In analogy with the LS-part, we now require that 
£> (NLS) (Q 2 ) pt be such as to absorb all the leading-/? parts of the difference (1/sf )[D(Q 2 )-V^ (Q 2 )] pt . In completely 
analogous way as before, we can show that this is equivalent to 

(-1)"( In" Ti + In"" 1 T x lnT 2 + • • • + In" T 2 ) (2) = (n = 1,2 . . .) , (B13) 

where X, = tjQ 2 e c / [i 2 , and coefficients c4n are defined in Eqs. (|52|l - (|54|) . These coefficients are known if the pertur- 
bative coefficients dj (|4"0|) are known. The (nonpower) expansion in Ak = -4fc(/ii 2 ) of the NLS-term is then 

s?p (NLS) (g2) = c W{ A l +AlA2 oC W +AlAa [^(2) +Ac P)] (B14) 

+ [^-^i^ 3 ]/3o(lnTilnT 2 ) (2) +0(^ 1 ^4,^2^3,.--)} • (B15) 

The last term in brackets has a coefficient oc (ln7i In '7 2 )r 2 ) which cannot be obtained on the basis of the perturbative 
coefficients dj (|4H)) . The perturbative part of this last term is zero. We know c^q if we know the NLO coefficient di 
of the perturbation expansion (|4"0|) of observable T>(Q 2 ); for the knowledge of we need, in addition, the knowledge 

(2) 

of <i 2 , and for c 22 the knowledge of e^. 

We now continue analogously one step further. In view of the normalization conditions (|B4|1 . it follows immediately 

1 

/V 

which is identical to the coefficient k^ (f57|) in the approach of Sec. IIVI We require that the third (N 2 LS) term 

S V V (K LS)(q2) in skeleton 

expansion (|B2j) satisfy the condition: 2?( N LS ^{Q 2 ) P t be such as to absorb all the leading- 
00 parts of the difference (l/sf)[V(Q 2 ) - V (LS \Q 2 ) - sf £> (NLS) (<3 2 )] P t- This then implies 

(- In T x - In T 2 - In %) (3) = c$ , (B17) 

where Cjj is given in Eq. (|60[) ; and similarly for higher terms (c 22 \ etc.). The (nonpower) expansion in Ak = ^(/Lt 2 ) 
of the N 2 LS-term is then 

s v V ^\Q 2 ) = s f{Al + A 2 1 A 2 0o^ 1 ) + O(A 2 A 3 ,A 1 Al...)} . (B18) 

We know the quantity if we know the coefficients d\ and d 2 in the perturbation expansion (|40p of observable 

V(Q 2 ); for the knowledge of we need, in addition, the knowledge of d%. 

Normalization conditions (|B4[) now imply that the coefficient of the N 3 LS-term in the skeleton expansion (|B1[) - 
dB2l is 
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which is identical to the coefficient k^ in the approach of Sec. IIVI The (nonpower) expansion in Ak = AkifJ 2 ) of 
the N 3 LS-term is then 

s ^2?( n3ls )(Q 2 ) = a fA$ + 0{A\A2,A\Az, ...). (B20) 



We know the quantity if we know the coefficients d%, d 2 and ds in the perturbation expansion (|40[) of observable 

Finally, we can combine the LS-term (|68[) . whose characteristic function is usually known, with all the beyond- 
the-LS terms written hitherto (|B14[) - (|B20[) which are known if d±, d 2 and d^ are known; since each of these terms 
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is RScl- independent, we can use in the most general case various (space-like) RScl's Q? = Q 2 exp(C 3 ) as Eq. (|A2|) 
(j = 2,3,4 for the NLS, N 3 LS and N 3 LS terms, respectively). This then results in 

V = p( LS >(Q 2 ) + 4 2) [A i m] 2 + {t i i ) A 1 (Ql)A 2 (Ql) + t^[A 1 (Qi)} 3 } 

+ {4 2 U 1 (g 2 M 3 (Q^) + 4 3) [^i(Qi)] 2 ^ 2 (Q§) + 4 4) [MQl)] 4 } 

+ { [A2{Ql)} 2 - A^QDA^Ql)} c^^anTj lnT 2 } (2) + 0{A\,A\A 2 , . . .) , (B21) 

where the coefficients t?l are precisely those given in Appendix [XI Eqs. (|A3[) - (|A14[) . Therefore, the evaluation method 
presented in the present Appendix, which is a representation of an assumed skeleton expansion (|B1|) - (|B2|) . reduces to 
the evaluation method presented in Sec. If VI when the following replacements are made: 

[A 1 {^)] kl [A 2 (n 2 )] k2 ■ ■ ■ [A s (n 2 )] k " i-> A kl +2k 2 -+sk 3 {t j2 ) ■ (B22) 

In the present method, the coefficient at the last term in brackets in expression (|B2f I) can be evaluated only if certain 
assumptions about the NLS characteristic function i^(ti, i 2 ) are made. For simplicity, we will make the factorization 
assumption 

F%(tiM) = w £ v {ti)w £ v {h) => (lnT 1 lnT 2 ) (2) = 1 (c^)' , (B23) 

where the last identity is obtained on the basis of identity (|Bf 3|) for n = 1. 

Similarly as the skeleton-motivated evaluation method (f70ll - (|7I|) . the skeleton evaluation method (|B21|) can be 
performed in principle at any chosen RScl's Qj and in any RSch of the class (|4Tj) . This method was denoted as 'vl ' 
in Ref. [12J. However, skeleton method (|B2f j) makes sense only if the skeleton expansion (|B2I) really exists. If the 
latter exists, it probably does so only in a specific ('skeleton') scheme [34], HH- In contrast, the skeleton-motivated 
evaluation method (|7D|) - ([7T|) does not rely on the existence of the skeleton expansion. 



APPENDIX C: LEADING-SKELETON CHARACTERISTIC FUNCTIONS IN THE SPACE-LIKE AND 

TIME-LIKE FORM 



In this Appendix we summarize the knowledge of the LS characteristic functions for the space-like observables 
V(Q 2 ). In the space- like formulation (|68|) . which involves the space- like coupling Ai, the characteristic function can 

be obtained from the knowledge of the leading- j3o coefficients Cnl - cf. Eqs. (|40|) and (|42|) . following the formalism of 
Neubert [H|. 

For example, in the case of the Bjorken polarized sum rule (BjPSR) db(Q 2 ), the leading-/3o coefficients were obtained 
in Ref. In MS RSch and at RScl /i 2 = Q 2 exp(C) (we use A = A throughout, i.e., C = C = -5/3) , they are 
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This implies that the (leading- (3q) Borel transform is 
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The renormalon poles are at u — ±1, ±2. The LS characteristic function appearing in (|68|) is then obtained by the 
general formula 
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where uq is any real number closer to the origin than the leading renormalon (—1 < uq < I). We can choose uq = 
and introduce a new integration variable r = — iu. The integral, with S(u) of Eq. (|C2j) . then reduces to 
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which can be calculated by the use of the Cauchy theorem in the complex r-plane: when r > 1, we close the path 
with a large semicircle in the upper half plane; when r < 1, in the lower half plane. This gives us the result 




which we already used in [12j. 

The LS characteristic function for the Adler function d,,(Q 2 ) was obtained in Ref. [3l{ . on the basis of the large- /3q 
expansion of the Borel transform of d v obtained in Refs. [5y, [57| 
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where C F = (N 2 - 1)/(2N C ) = 4/3. 

The semihadronic r decay ratio r T is a time-like observable. The LS term of r T can be obtained from the LS-term 
of the Adler function on the basis of the relation 



r T (AS = 0,m q = 0) = — 
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This implies for the LS term of r T (AS = 0, m q = 0) 

r T (AS = Q,m q = 0)< LS ) = / ^ F r M (t) Sli^m 2 .) , (C9) 

Jo t 

where 21 1 is the time-like coupling appearing in Eqs. (|T])- (|10p . and superscript A4 in the characteristic function means 
that it is Minkowskian (time-like). The latter was obtained in Ref. [32j| 10 
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Here, PolyLog„ is the polylogarithm function of n'th order (using notation of |13j). 

The LS part of any space-like observable T>(Q 2 ) can be written in two equivalent forms - the form involving the 
space-like coupling A\, Eq. (|68p. and the form involving the time-like 2li of Eqs. ([7|- (fT0l) 
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Jo 1 Jo 1 
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where the superscript M. stands for the "Minkowskian" (time-like) formulation, and the two characteristic functions 
arc related via relations 
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We use a different normalization, so an additional factor of t/4 appears, in comparison to [32(1 . 
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Identity (|C14[) is a direct consequence of the definition of J--p there. On the other hand, relation ()C12|) is a direct 
consequence of identity (|(J13|) and of the following identity in the complex cr-plane (where a — k 2 is square of a 
four- vector) : 

/ — A^K 2 = -ae c ) [T v (<i/Q 2 ) - T v (0)] = , (C15) 
JC1+C2 a 

where function Tx> is defined by identity (|C14|) , and the path C\ + C 2 is depicted in Fig. ([8]). In the cx-plane, the 



a- plane 



c 2 



FIG. 8: The path of integration of integral JC15ft in the complex cr-plane: a > semiaxis is the cut of *4i(— 0"e c ) factor, and a < is the 
cut of the [J r T>(o'/Q 2 ) — ^dCO)] factor in the integral. 

only singularities of the integrand in Eq. (|C15[) are the cut of A\(— a) along the positive semiaxis, and the cut of 
\Txs{p I Q 2 ) ~ J--d{0)] along the negative semiaxis. Identity (|C15[) thus follows from the Cauchy theorem. 

When applying relation (|C13|) to the characteristic function (|C5j) of BjPSR, we obtain for the time-like characteristic 
function of that observable 



Frit) = * 



10 122 2/ 21 



(C16) 



which agrees with the corresponding expression in Ref. [58| after identifying in their Eq. (4.57): ^(e, N = 1) = 
(-3/2)F 6 M (t), and e ee t = p?jQ 2 (Ref. uses apparently C = 0). 

APPENDIX D: EXPLICIT EXPRESSIONS FOR VARIOUS COEFFICIENTS 

Expansion @ is solution of the perturbative RGE equation (jTJ) . If the conventional ( "MS" ) reference scale A fLU Il5| 

is adopted, and RGE |T]) is iteratively solved for large Q 2 /A [ln(Q 2 /A ) S> 1] in an arbitrary RSch (/3 2 ,/33, ■ • •)) this 
results in expansion ^ with coefficients Km given in Eqs. ([3]) for k < 3, and for k = 4, 5, 6 given below (notations: 
c,- = &/#,): 
For k — 4: 

*« = -Kl) ('-!)■ K —(l)(- 2+3 ? 



5(4 



^ = 5(3). *«3 = -(£)- (Dl) 



For fc 



i^so = i (7c* - 18c 2 c 2 + 10c 2 - cic 3 + 2c 4 ) , 

= 4 (4c? - 3cic 2 - 2c 3 ) , K 52 = (cf-4c 2 c 2 ) , 
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For k — 6: 



K 60 = — — g (l7c? - 18c 3 c 2 - C1C2 - 23cjc 3 + 24c 2 c 3 - 2cic 4 + 3c 5 ) , 



1 

I 

5 _i_ 07^3, 



^61 = (-He? + 72cf c 2 - 50cicl - 7c?c 3 - 10cic 4 ) 
1 



K 62 = -rg (~23cf + 27c?c 2 + 10c? c 3 ) 



1 



77cf 



^64 - tt^, ^65 = -^ . (D3) 



In practical calculations, we use: (a) at loop-level=3: C3 = C4 = C5 = and we include in expansion J2J terms if/^ up 
to fc max = 5; (b) at loop-level=4: C4 = C5 = and we include terms up to fc max = 6. 



The perturbation coefficients dj (j = 1,2) of the perturbation expansion for the massless Adler function d v (Q 2 ), 
cf. Eq. (J40J), in MS RSch and at RScl \i 2 — Q 2 , are known exactly, Refs. [H,[33], respectively 

d[ AdL) = — + 0.691772/3 , d { 2 Adh) = -27.849 + 8.22612/3 + 3.10345/3q • (D4) 

The N 3 LO coefficient cfe, in the aforementioned RSch and RScl, was obtained in an approximate form in Ref. [38| 
[(Eqs. (20) and (12) in 38])]: 

4 Adl ° = 46.1992- 131. 04/3 + 49. 5237/^ + 2. 18004/3^ , (D5) 

where the coefficients at (3q and at /3g are known exactly ([{36|, H3], [IB]), and the other two coefficients were estimated 
in Ref. [38| by using the methods of the principle of minimal sensitivity (PMS) (29j, and of the effective charge (ECH) 

urn. 

The light-by- light contributions are not included in the coefficients (ID4I and (|D5[) . They have a different topology 
of diagrams and should probably be resummed separately (cf. Ref. jll|), and they appear for the first time at ~ a 3 
[37l ] . They are proportional to the sum of the charges J2Qf- This sum is zero in the case n/ = 3 considered here. 

Coefficients d\ and di for BjPSR db{Q 2 ), in the aforementioned RSch and RScl, were obtained in Ref. [4l[ and are 

4 Bj '° = -— + 2/3 , 4 Bj ° = -35.7644+ 10. 5048/? + 6.38889/3^ . (D6) 

In the coefficient d^'\ only the leading-n/ part (cx n 3 ) is known exactly [hj (•o- the leading-/3 part, oc /3q). On 

this basis, the authors of Ref. [12] obtained estimates of 4 as a polynomial in (3q by using naive nonabelianization 
(NNA) [43]: n/ 1— > — 6/3o. Several relations between BjPSR, Bjorken unpolarized sum rule, and Gross-Llewellyn Smith 
sum rule were found out and investigated in Refs. [6l|. 

APPENDIX E: MASSLESS PART OF THE STRANGELESS TAU DECAY RATIO 

In this Appendix we extract the measured value of the massless part of the QCD-canonical strangeless ratio 
ry(AS = 0, m q = 0) for the semihadronic decay, on the basis of the results of the final ALEPH data analysis [j| HI]. 11 
This quantity is related to the ALEPH-measured [1, |j| (V+A)-decay ratio 

m*s=o) . - RA as^) ( ed 

F(t -> v T e Veil)) 



(1-B e - 



-Rr(AS^O) = 3.482 ±0.014 . (E2) 



For an extraction of r T (AS = 0, m q = 0) based on the older set of measured results 0, see for example Ref. l62l . 
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These values were obtained in Ref. @ from the measured leptonic branching ratio B e = B(t — > v T e v e ) = (17.810 ± 
0.039)% (ALEPH, [|), from = B(t~ -> v^V^) = (17.332 ± 0.049)% (world average, @), and from the 
strangeness-changing branching ratio B$ = (2.85 ± 0.11)% (ALEPH, fH). The relation between the canonic massless 
quantity r T (AS = 0, m q = 0) and quantity (|E1|) - (|E2|) is 

r T (AS = 0,to,j = 0) = r T (AS = 0,m q ) - 5r T {AS = 0,m M ^ 0) (E3) 

R T (AS = 0) 



3|Kd| 2 (l + few) 



(1 + <W) - Sr T (AS = 0,m Uld ^ 0) . (E4) 



Here, r T (AS = 0,7n q = 0) is QCD-canonical, i.e., its pQCD expansion is r T (AS — Q,m q = 0) pt = a + 0(a 2 ); the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix element \V u d\ has the value largely dominated by + — > + nuclear beta 
decays [63j 

\V ud \ =0.9738 ±0.0003 , (E5) 

the electroweak (EW) correction parameter is 1 + 5ew = 1-0198 ± 0.0006 [1, Q; the residual EW correction parameter 
is S' EW = 0.0010 [HI; the (V+A)-channel corrections Sr T (AS — 0,m u _d ^ 0) due to the nonzero quark masses are 

fijl. |65| the sum of the D = 2,4, 6, and 8-dimensional corrections (S^ v and their value is @ either 5r T (AS = 

0, m U: d 7^ 0) = (— 5.2±1.7) x 10~ 3 if the gluon condensate contribution is included, and (— 5.0±1.7) x 10~ 3 if the gluon 
condensate contribution is not included (using for the gluon condensate the ALEPH-value (aGG) = (— 0.5±0.3) x 10~ 2 

Bi). 

Inserting all the aforementioned values in relation (|E4[) and taking into account the value (|E2|) . we extract the 
experimental prediction for r T (AS = 0, m q = 0) based on the most recent ALEPH data 

rv(AS = 0, m g = 0) cxp . = 0.204 ± 0.005 , (E6) 

where the uncertainties have been added in quadrature. The uncertainty in Eq. (|E6p is dominated by the experimental 
uncertainty 5R T = ±0.014 (|E2[) . The value (|E6|) remains unaffected up to the displayed digits when we either include 
or exclude from the above quantity the gluon condensate contribution. 
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